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Abstract 



We perform an average analysis of the Grassmann condition number ^{A) for the 
homogeneous convex fcasibihty problem 3.T G C \ : Ax = 0, where C C M" may be any 
regular cone. This in particular includes the cases of linear programming, second-order 
programming, and semidefinite programming. We thus give the first average analysis of 
convex programming, which is not restricted to linear programming. The Grassmann 

1-^^ condition number is a geometric version of Renegar's condition number, which we have 

introduced recently in [arXiv:1105.4049vl]. In this work we use techniques from spherical 
pH convex geometry and differential geometry. Among other things, we will show that if the 

I— I entries of A G jj™^" are chosen i.i.d. standard normal, then for any regular cone C we 

have E[ln'^(A)] < 1.51n(n) + 1.5. 

> 

CD AMS subject classifications: 90C25, 90C22, 90C31, 52A22, 52A55, 60D05 

^ Key words: convex programming, perturbation, condition number, average analysis, spher- 

Cn| ically convex sets 

-H 1 Introduction 



Convex programming is an efficient tool in modern applied mathematics. In fact, a commonly 
accepted technique in current scientific computing is to "convexify" supposedly hard problems, 
solve the relaxed convex problem, and then hope that the result is close to a solution of the 
original problem. To quote from [9, §1.3.2]: "With only a bit of exaggeration, we can say that, 
if you formulate a practical problem as a convex optimization problem, then you have solved 
the original problem." 

But what is the complexity of convex programming? To specify this question further, 
we ask for the number of arithmetic operations, or the number of iterations of an interior- 
point method. Steve Smale suggested in [49] to use the concepts of condition numbers and 
probabilistic analysis in a 2-part scheme for the analysis of numerical algorithms: 1. Establish 
a bound for the running time, which is polynomial in the size of the input and (the logarithm 
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of) a certain condition number of the input. 2. Analyze the condition number of a random 
input in form of tail estimates. 

The first step of this scheme, i.e., the analysis of the role of condition numbers in convex 
programming, was initialized by Jim Renegar in [42, 43, 44], and is an active area of research, 
cf. [55, 56, 23, 26, 25, 39, 41, 17, 18, 21, 40, 24, 54, 6]. In these references the role of condition 
numbers is analyzed for linear and nonlinear convex programming, for exact arithmetic and 
for finite-precision arithmetic, for ellipsoid and for interior-point methods, etc. 

Yet, the second step of Smale's scheme, i.e., the probabilistic analysis of the condition 
number, was until now severly restricted to the linear programming case. See the survey arti- 
cle [10] and the references given therein for more details on probabilistic analyses of condition 
numbers for linear programming. 

We will give in this paper the first average analysis of a condition number for the general 
homogeneous convex feasibility problem. This includes the special cases of linear program- 
ming, second-order programming, and notably also the semidefinite programming case. More 
precisely, we consider the following problem: 

Let C C be a regular cone, i.e., C is a closed convex cone with nonempty interior that 
does not contain a nontrivial linear subspace. The dual cone^ of C is defined as C := {z e 
R"" I Vx G C : z'^x < 0}. We call C self-dual if C = — C. The homogeneous convex feasibility 
problem is to decide for a given matrix A G M"*^", 1 < m < n, the alternative^ 

3x G M" \ s.t. Ax = 0, xeC , (P) 
ByGM^'VO s.t. A^yeC. (D) 

This problem reduces to the linear feasibility problem if C = , it reduces to the second-order 
feasibility problem if C = £"1 x . . . x where £."■:= {x eW \ Xn > xj + ■■■ + ... , xl_-^} 
denotes an n-dimensional Lorentz cone, and it reduces to the semidefinite feasibility problem 
if C = Sym^ = {M G Sym^ | M ^ 0} is the cone of positive semidefinite matrices, where 
Sym'^ = {M G W'"'"' \ = M}. 

The condition number for which we will provide an average analysis, is the Grassmann 
condition number, that we have introduced in [3, 4], cf. also [6]. In the following paragraph we 
will recall the necessary definitions so that we can state the main theorem of this paper. See [4] 
for a more extensive description of the Grassmann condition and its relation to Renegar's 
condition number. We denote by GTn,m the set of m-dimensional linear subspaces in M". It is 
a well-known fact that Gr^j m is a smooth Riemannian manifold, called Grassmann manifold. 

Definition 1.1. Let C C M" be a regular cone, and let 1 < m < n — 1. We define the sets of 
dual feasible and primal feasible subspaces via 

Vm{C) := {W G Gr„,^ \WnC^O}, Vm{C) := {W G Gr„,„ \W^nC^O}. 

Furthermore, we define the set of ill-posed subspaces via 

'^Some authors call C the polar cone, and define the dual cone as —C. 

^In fact, (P) and (D) are only weak alternatives, as it may happen that both (P) and (D) are satisfiable. 
But the Lebesgue measure of the set of these ill-posed inputs in R""*" is zero. 
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It is known (cf. [4]) that V^iC) and Vm{C) are compact subsets of Gr„^m, and Vm{C) U 
'Pm{C) = GTn,m- Moieover the boundaries of T>^{C) and VmiP) coincide with E^(C). Fur- 
thermore, 

E„(C) = {W e GTn,m I n C ^ and W n int(C) = 0} . 

In other words, the set of ill-posed subspaces consists of those subspaces, which touch the 
cone C at the boundary. The following definition is a restatement of [4, Def. 1.2]. We denote 
by dg the geodesic distance in Gr„_j„ (see [4] for more details on this). 

Definition 1.2. The Grassmann condition number is defined via 

% : G.„,„ ^ [1, oo] , ,^,^^^^^^^C)) • 

where dg{W,^m{C)) := mi{dg{W,W') \ W G S„,(C)}. 

To simplify the notation we write ^{W) := ^c{W)- Furthermore, for A G M*"^" with 
i]<i{A) = m we write ^{A) := ^{im.A'^). Before we state the main result of this paper, we 
recall from [4] the following fundamental relation between the Grassmann condition number 
"^{A) and Renegar's condition number Tl{A) 

"^{A) < n{A) < KiA) ■ '^{A) , (1) 

where k{A) denotes the usual matrix condition number, i.e., the ratio between the largest 
and the smallest singular value of A. The Grassmann condition number can be interpreted 
as a coordinate-free version of Renegar's condition number, and the dependence of Renegar's 
condition number on the representation of the subspace is conveniently described by the matrix 
condition number k{A). See [4] for more details. 

The inequalities (1) show that probabilistic analyses of ^ easily imply corresponding results 
about TZ, as the probabilistic behavior of the usual matrix condition number is a well-studied 
subject, cf. [20, 16, 13]. For example, in [16] it was shown that for Gaussian matrices A G M"*^" 
one has 

71 

E\lnK(A)] < ln( , r) + 2.3 . (2) 

^1 + \n — 

The following two theorems are the main results of this paper. 

Theorem 1.3. Let C C 6e a regular cone with n>3. If A e M"*^" is a Gaussian random 
matrix, then we have 

1 3 

Prob[^(y4) > f] < 6 • y/m{n - m) ■ - , iftyn^^, (3) 
E[ln^(^)] < 1.5 • ln(n) + 1.5 . (4) 

In fact, these bounds hold for any probability distribution on M"^^", which induces the 
uniform distributions on the Grassmann manifolds Gr„^„_m and Grn,m via the maps A i— t- ker A 
and A i— )• imA^. We have reduced our formulation to the case of Gaussian random matrices 
for simplicity. 

Assuming certain conjectures on the growth of the intrinsic volumes of special cones C 
(see Section 2.3) we can improve the bounds as follows. 
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Theorem 1.4. Let C C M" 6e one of the following cones: M.^, C, £"1 x . . . x Sym^. 
If A E M"*^", with m > 8, is a Gaussian random matrix, then we have 

Prob[^(A)>t] < 20-v{C) ■ Vm-j , ift>m, (5) 
E[ln^(^)] < ln(m) +max{ln(t;(C)),0} +3 , (6) 

where the excess over the Lorentz cone v(C) (cf. (28) for the definition) is hounded as shown 
in the following table 
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To illustrate the significance of Theorem 1.3, we will use a result from [54] to state a 
complexity result about general convex programming. 

Corollary 1.5. Let C he a self-scaled cone with a self-scaled harrier function. Then there exists 
an interior point algorithm, that solves the general homogeneous convex feasihility prohlem 

for Gaussian random inputs in expected O^^fuc ■ (ln(z/c') + ln(n))) number of interior-point 
iteraMons. Here, vc denotes the com,plexity parameter of the harrier function for the reference 
cone C . For linear programming, second-order programming, and semidefinite programming, 
we have an expected number of interior-point iterations of 0{^\fn- ln(n)). 

Proof. In [54] the authors describe an interior-point algorithm that solves the general homo- 
geneous convex feasibility problem, for C a self-scaled cone with a self-scaled barrier func- 
tion, in 0{y/iyc • ln(/^c' ' ^(^))) interior-point iterations. From (1) we have \n{TZ{A)) < 
ln{K{A)) + ln(^(A)). Using the estimates (2) and (4) wc thus get E[ln(7^(^))] = 0(hi(n)). 
As for the complexity parameter uq, the typical barrier functions for (LP), (SOCP), (SDP), 
yield 

C = Rl : uc = n , C = C"' x . . . x C^- : uc = 2r , C = Sym^ : uc = k . 
In particular, for these cases we have vc ^n. □ 

Additionally, we mention that in [54] it is also shown that the condition number of the 
system of equations, that is solved in each interior-point iteration, is bounded by a factor 
of 7^(^)^. One can thus also use our results to get information about the expected cost of 
each iteration in the above-mentioned algorithm. We leave this to the interested reader. 

This paper overall consists of two parts. The first part, which consists of the Sections 2- 
3, is elementary and easily accessible. In Section 2.1 we review some facts about spherical 
intrinsic volumes and in Section 2.3 we explain how these may be estimated. In Section 3 we 
use the fact that Gr^^^ is a Riemannian manifold, which allows to define the tube TiT^m-, oi) 
of radius a around with respect to the geodesic distance dg, cf. (33). We will state in 
Theorem 3.1 an estimate of the volume of T{Tirn-i aiid we shall derive proofs for Theorem 1.3 
and Theorem 1.4 from this estimate. 

The second part of the paper is devoted to the proof of Theorem 3.1. In Section 4 we will 
define a certain generalization chy ((^, t) of the characteristic polynomial of a self-adjoint linear 
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operator : V ^ V oi a euclidean vector space V, which depends on a hnear subspace Y CV. 
We call chylfjt) the twisted characteristic polynomial and we will compute its expectation 
when the subspace Y (of a fixed dimension) is chosen uniformly at random. Section 4 only 
requires some elementary linear algebra, and although it seems to lie outside the general 
framework of this paper, it is a central step towards the proof of Theorem 3.1. In Section 5 
we will provide some preliminaries from Riemannian geometry, and finally, in Section 6 we 
will complete the proof of Theorem 3.1. 

Prom a high level point of view, one may interpret Theorem 3.1 as an extension of Weyl's 
(spherical) tube formula, cf. [57], to the Grassmann manifold setting. Although such extensions 
in the sense of Theorem 3.1 were known before, cp. [28, 29], these have the serious drawback 
that they only hold for radii a below a certain threshold qq, which depends on the cone C. In 
fact, for the interesting cones used in convex programming, this threshold is ao = 0, so that the 
previously existing results are useless for applications in convex programming, cf. Remarks 2.1 
and 3.2. Theorem 3.1 does not suffer from this restriction and may thus be interpreted as 
a less precise, but more robust version of Weyl's tube formula in the Grassmann manifold 
setting, like it was achieved for the spherical setting in [14]. 

A key idea in the proof of Theorem 3.1 is to consider cones with smooth boundaries. More 
precisely, we will assume that M := dC H S^~^ is a smooth hypersurface of S"""^ with strictly 
positive curvature. Then it turns out that the corresponding set ^^(C) of ill-posed inputs 
is an embedding of the Grassmann bundle Gr(M, m — 1) in the Grassmann manifold Gr„^^. 
Exploiting the induced bundle structure on is the crucial idea for the proof of Theorem 3.1. 

2 Spherical convex geometry 

It is a simple, yet essential observation that the analysis of the question whether a subspace 
hits a convex cone nontrivially, may be transferred from euclidean space to spherical space by 
intersecting both the subspace and the cone with the unit sphere. The above question then 
translates into the question of whether a subsphere of the unit sphere intersects a spherical 
convex set, or not. The analysis of the homogeneous convex feasibility problem thus naturally 
finds its place in the domain of spherical convex geometry. 

While euclidean convex geometry is a classical and extensively studied subject, the situa- 
tion for spherical convex geometry is quite different. The theory here is much less established, 
while a number of results, which are difficult to find in the literature, seem to be folklore 
among the experts in convex geometry. We find it thus appropriate to recall in this section a 
number of elementary facts and definitions to provide a preferably self-contained presentation. 
In our presentation we rely on the theses [28] (cf. also [29]) and [3], and the survey article [27]. 
We will put a special emphasis on intrinsic volumes. 

The first subsection is mainly devoted to the elementary definitions and the corresponding 
notations. In the second subsection we will give the definition of spherical intrinsic vohimcs, 
and we will present formulas for the intrinsic volumes of the cones 

In the third subsection we will describe a general way to estimate the intrinsic volumes (of 
self-dual cones). These estimates will be used in the proof of Theorem 1.4, which we will give 
in Section 3. 
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2.1 Some general facts 

We denote the spherical distance function by d: S'^~^ x S'^~^ M, d{p,q) = aTCCOs{{p,q)). 
For a subset K C S^~^ of the unit sphere we denote the cone generated by this set by 
cone(i^) := {A • p I A > , p € K}, and we additionally set cone(0) := {0}. A subset 
K C S"'~^ is called (spherical) convex iff for all p,q ^ K with q / ±p the (unique geodesic) 
arc between p and q is contained in K. This is equivalent to the condition that the set cone(if) 
is a convex cone. We denote the family of closed spherical convex sets by /C(S'"-^). Note that 
the empty set as well as the whole sphere S""^"*^ both lie in K-^S'"^^^). 

For M C 5"^-^ closed and a > we denote the tube of radius a around M by T{M, a) := 
{p G S'^"^ I 3g G M : d{p,q) < a}. Using the notion of tubes we may define the Hausdorff 
distance dn on ]C{S"'''^) in the following way 

d^{Ki,K2) :=max{min{a > \ K2 C T{Ki,a)} , mm{P > | i^i C r{K2, ^)}} , 

for Ki,K2 G IC{S^~^). It is well-known (cf. for example [38, §1-2]) that the map d^ is a metric 
on ]C{S"'~^), which turns ]C{S"'~^) into a compact metric space. 

The duality map which sends a closed convex cone C onto its dual cone C naturally defines 
an involution on IC{S^~^), which we also denote by /C(S'"'~^) — t- IC{S^~^), and which is given 
by K := (7n S'^-i = {p ^ S''-^ \ d{p,K) > f }, where C := cone{K). It can be shown (cf. for 
example [28, Hilfssatz 2.2] or [3, Prop. 3.2.4]) that for -fCi, K2 G /C(S"'-i) with duiKi,K2) < f 
we have du{Ki, K2) = dii{Ki, K2). In particular, the duality map is a local isometry. 

Remark 2.1. One important difference between euclidean and spherical convex geometry are 
the convexity properties of tubes. In the euclidean case the tubes around a convex set are 
again convex. In the spherical case it is rarely true that the tube around a convex set is again 
convex, unless the tube is the whole sphere. For example, it is elementary to show (cf. [3, 
§3.1.2]) that ii K G }C{S"'~^) contains two antipodal points, or if the boundary of K contains 
a geodesic segment, then T{K,a) G ]C{S"'~^) for some a > implies T{K,a) = S""^. It 
follows that tubes around (proper) subspheres are either the whole sphere or not convex. The 
same statement holds for spherical convex sets K, whose corresponding cone C = cone{K) 
has a nontrivial face, i.e., there exists a supporting hyperplane H C M", H Ci int(C) = 0, such 
that the (linear hull of the) face F := H d C has dimension dimlin(F) > 2. Hence, among 
our standard ensemble of cones R" , £"1 x . . . X Sym^, only the cones n G N, have 
convex tubes different from the whole sphere. 

Another difference between euclidean and spherical convex geometry is that while the set 
of all convex bodies, i.e., the set of all compact euclidean convex sets, is connected, this is not 
true for /C(5'"~^) as we will see next. We use the notation 15^^(5*"^^) := {0}, and 

5^(5"-^) := {S C 5"-^ I 5 is a fc-dim. subsphere} , A: = 0, 1, . . . , n - 1 , 
JCiS""-^) := {K G /C(5"-^) \K^$ andKis not a subsphere} . 

Note that there is a canonical bijection between 5*^(5""^) and Gr^^k+i given by Gr^ ^^i — >■ 
5fc(5"-i), w^Wn and its inverse cS'=(5"^^) ^ Gr„,fc+i, S ^ lm{S). Note also that 

S = S^ :=W^n 5"-^ for S a subsphere of 5"-^, with W := lm{S). 

As the non-subspheres are central objects for our study, we call elements in }C^{S"'~^) caps. 
This naming is different from other works, where "cap" may denote a spherical ball (cf. for 
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example [11]). We will denote spherical balls by the term circular caps. It is easily checked 
that an element K G /C(S'"~^) is a cap iS 3p £ K : —p ^ K. 

The set of spherical convex sets K,{S'^~^) decomposes into caps and subspheres, i.e., the 
decomposition of /C(5""^) into its connected components is given by 

n-l 

/C(5"-^) = /C"(5"-^) U (J 5*^(5"-^) . (7) 

k=-l 

The connectedness of the components in this decomposition is checked easily. The fact that 
elements from separate components cannot be connected via a path in K{S'^~^) is seen in the 
following way. If Wi , W2 are subspaces of M" of different dimensions, then counting dimensions 
shows that Wi n 7^ or fl W2 / 0. This implies that the corresponding subspheres 
Si := Wi n i = l,2, have Hausdorff distance d^{Si,S2) > f . In fact, as T{S, f ) = 5""^ 

for any (nonempty) subsphere S of we have dii{Si,S2) = f for subspheres of different 

dimensions. For a cap K G /C'^(5'"~^) and a subsphere S we have the well-known theorem 
of alternatives, also known as Farkas- Lemma, which says that K Ci S ox K Ci S-^ / 0. 
This implies dji{K,S) > |. The decomposition in (7) is thus the decomposition ]C{S"'~^) in 
connected components. 

The set of caps /C'^(S'"^^) should be seen as the essential part of /C(5"^^) containing a 
variety of sets with diverse properties. We need to specify subfamilies of /C(5"^^) with which 
we can work in a unified way. We begin with the set of polyhedral convex sets: A spherical 
convex set K G ]C{S'^~^) is called polyhedral if cone(if) is the intersection of finitely many 
half-spaces: 

^P(^n-i-) ._ 1^ g /C(5"-^) I cone(ii:) = Hin...nHk, Hi is a closed half-space in M"}. 

Equivalently (cf. for example [45, Sec. 19]), a set K G /C(5"'^) is polyhedral iff it is the 
(spherical) convex hull of finitely many points in 5""^. Our standard example for a polyhedral 
convex set is the intersection M" n S^~^ of the positive orthant with the unit sphere. By 
definition, the set of polyhedral convex sets contains the subspheres of S*"^^. Furthermore, 
the set of polyhedral convex sets lies dense in /C(>S'"'~^) with respect to the Hausdorff metric. 
This is seen by an easy adaption (cf. [28, Hilfssatz 2.5] or [3, Prop. 3.3.4]) of the proof for the 
euclidean statement (cf. [47, §2.4]). 

Another important subfamily of spherical convex sets is given by the set of regular caps: A 
spherical convex set K G /C(5"~^) is called regular if both K and K have nonempty interior: 

/C'-(5"-i) := {K G /C(5"-i) I int(i4:) 7^ and mt{K) 7^ 0} . 

Note that subspheres are not regular, i.e., the set of regular convex sets is a subset of the set 
]C'^{S^~^) of caps, and we may thus speak of regular caps. In the set /C^(S'"~^) of regular caps 
we define for n > 3 the subclass of smooth caps via 

dK is a smooth hypersurface in S"'~^ \ 
with nowhere vanishing Gaussian curvature J 

Here, the Gaussian curvature is the determinant of the Weingarten map (see for example [53] 
for an elementary introduction of these notions; see [19, Ch. 6] for a more thorough treatment). 
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For hypersurfaces M of S'" with a unit normal vector field v. M ^ S'^ the Weingarten 
map at p G M is given by 

Wp : TpM ^ TpM , WpiQ = -Dpv{Q , (8) 

where DpU denotes the derivative of u at p (for hypersurfaces of M" this is elementary (cf. for 
example [53, Ch. 9]); for hypersurfaces of S"'~^ this is also true (cf. [3, Sec. 4.1.1])). It 
can be shown that the Weingarten map is self-adjoint, and we will denote its eigenvalues by 
Ki{p),. . . ,Kn-2{p)- These are called the principal curvatures of M at p. Furthermore, we will 
denote hy ak{p), < k < n — 2, the fcth elementary symmetric function in the eigenvalues 
of Wp. When M is the smooth boundary of a cap M = dK with int(ii') / 0, then we 
will always assume that u{p) denotes the normal direction pointing inwards the cap K. This 
has the effect that the Weingarten map is positive definite, i.e., M has everywhere positive 
curvature (cf. [47, §2.5] for the corresponding euclidean statement). 

Example 2.2. Let M = dK C S'^~^ be the boundary of a circular cap K := B{z, (3) = {p e 
S"^-*- I d{z,p) < f3} of radius /? G (0,7r). Then the principal curvatures Ki{p) , . . . , Kn-2ip) 
at p G M are given hy Ki{p) = ... = K„_2(p) = cot(/3). In particular, we have ai{p) = 
i^f) . cot(/3)\ 

For the euclidean case, Minkowski (cf. [7, §6]) has shown that the set of smooth convex 
bodies lies dense in the set of convex bodies. Using this result, it is elementary to show (see 
for example [3, Prop. 4.1.10]) that the set of smooth caps /C™(S'"~^) lies dense in the set of 
all caps /C'^(5"-i). 

We finish this section with an important property of smooth caps, that we will need in the 
proof of the main result in Section 6. 

Proposition 2.3. Let K G /C™(5"-i) and C := cone(i^). IfW e S„(C), then Wr\K = {p} 
for some p G dK. 

Proof. As G Sm(C), we have n C = n 5C 7^ {0}. It follows that there exists 
p ^W^K. To prove that p is the only element 'm.Wf^K, assume that there exists WHK, 
p ^ q. As is a regular cap, wc have p 7^ — g, so that there exists a unique geodesic arc 
between p and q. By convexity of W and K, this arc lies in fl i^, and thus in M, the 
boundary of K. But this implies that along this arc M has zero Gaussian curvature, which 
contradicts the assumption K G K^'^{S'^~^). □ 



2.2 Intrinsic volumes 

Before we give the definition of spherical intrinsic volumes we introduce the following notation 
for the volume of the unit sphere and for the volume of tubes around subspheres 

0„_, := vol„_,(S"-') = (9) 

On-l^k{oi) ■■= voln-iT(5,a) , 

where S G 5*^(5"^^), and < a < |. An elementary computation shows (see for example [12, 
Lem. 20.5] or [3, Prop. 4.1.18]) that the volume of the a-tube T{S,a) is given by 

0„_i,fe(a) = Ok ■ On-2-k ■ r cos{pf ■ sin(p)"-2-*^ dp . (10) 

Jo 
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In particular, the volume of a circular cap B{z, f3) of radius /3 G [0, tt] is given by 

vol„_i B{z, (3) = On-2 ■ / sin(p)"-2 dp . (11) 

Jo 

The following proposition may be used for the definition of the spherical intrinsic vol- 
umes (cf. [27]); it may be interpreted as a spherical version of the Steiner polynomial (cf. for 
example [47]). 

Proposition 2.4. For K G /C(S'"~^) \ {0} and < a < | the volume of the a-tube around K 
is given by 

n-2 

vol„_i r(K, a) = vol„_i(i^) + VjiK) ■ On-ij(a) , (12) 

3=0 

for some uniquely determined continuous functions Vj : ]C{S'^^^^) — t-M, < j < n — 2. 
Proof. See (one of) the references given in [27]: [31, 2, 46, 36, 28]. □ 

Definition 2.5. For — l<j<n — 1 the j'th (spherical) intrinsic volume is the function 
Vj : /C(5"~^) M, defined as follows. For < j < n - 2 it is defined via 

J the quantity defined in (12) if / , 
'^^^''^'■=[0 ifK = 0. ^''^ 

Furthermore, for j G {—1, n — 1}, the jth intrinsic volume of K is defined via 

f if 7f ^ (?) 

Vn-,{K):=l'^--^ ifif = 0' V.^{K)■.= Vn.,{K). (14) 

The definition of V-i{K) in (14) may seem artificial at first sight, but in fact it is natural, as 
win become clear in (15), (16), (17) below. We will use the notation Vj{C) := Vj{C D S""-^), 
if C C M" is a closed convex cone, and refer to Vj{C) as the jth intrinsic volume of the 
cone C. Furthermore, we use the notation V{K) for the (n+ l)-tuplc of the intrinsic volumes, 
i.e. V{K) := (F_i(i^), Fo(i^), • • • , K-i(^)), and similarly for V{C). As for the notation, it 
turns out that usually the formulas become nicer when using the shifted index Vj-i{C) with 
< j < n. However, we refrain from reindexing the functionals Vj to simplify the notation, 
as the given definition is by now established (cf. [28, 27, 48, 3]), and a redefinition could lead 
to major confusion. 

The following proposition provides a simple and illustrative description of the spherical 
intrinsic volumes in the special case of polyhedral cones. Additionally, a couple of essential 
properties of the intrinsic volumes are immediate from this description. 

Proposition 2.6. Let C C R" be a polyhedral cone and let He '■ — >■ C denote the projection 
map onto C. Furthermore, let dc denote the function 

: ^ {0, 1, 2, . . . , n} , dim(face(nc(x))) , 

where face(a;) denotes the face of C of which x lies in its relative interior. Then the (j — l)th 
intrinsic volume of C, < j < n, is given by 

lS_i(C) = Prob [dcip) = 3] = Prob [dcix) = j] , (15) 
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where p E S'^ ^ is drawn uniformly at random and x G is drawn at random according to 
the normal distribution A/'(0, /„). 

Proof. This was shown in [3, Prop. 4.4.6]. The statement is essentially a reinterpretation of a 
well-known formula for the spherical intrinsic volumes of polyhedral cones, cf. for example [27, 
Eq. (11)]. □ 

Remark 2.7. 1. Prom Proposition 2.6 it is immediate that the intrinsic volumes satisfy 

n 

yj = 0,...,n:Vj.,{K)>0, J2^j_,{K) = 1 (16) 

j=0 

for all K G /C(S'"~^). The vector V{K) may thus be interpreted as a probability distri- 
bution on the set {0, 1, . . . , n}. 

2. Another consequence of Proposition 2.6 is the following nice duality property of intrinsic 

volumes 

F,_i((7) = K_,_i(C), (17) 

where < j < n. This relation may be deduced (cf. [3, Prop. 4.4.10]) from the duality 
relation between the face lattices of C and of C, respectively. 

Remark 2.8. It is known that the Euler characteristic % on ]C{S"'^^) has the representation 

n 

j=0 
j odd 

(See [28, Sec. 4.3] or [48, Thm. 6.5.5] for proofs of this fact.) Since xi^) = 1 ^^r all caps 
K G /C'^(5"~^), it follows that the intrinsic volumes of caps are bounded by Vj-i{K) < ^. 

Another simple consequence of Proposition 2.6 is a calculation rule for direct products of 
cones. The euclidean analog of this statement is well-known (cf. for example [34, Sec. 9.7]); 
the spherical statement below was, to our best knowledge, for the first time given in [3]. 

Corollary 2.9. Let C\,C2 he closed convex cones. Then the intrinsic volumes of the direct 
product Ci X C2 are given by 

ViCi X C2) = V{Ci) * V{C2) , (18) 
where * denotes the convolution operator, i.e., Vfe_i(Ci x C2) = Yl ^-i(C'i) • Vj-i{C2). 

i+j=k 

Proof (Sketch). We sketch the main ideas of the proof and refer for a more elaborated proof 
to [3, App. B.l]. By continuity of the intrinsic volumes, it suffices to consider the case where 
both Ci and C2 are polyhedral cones. The faces of Ci x C2 are given by direct products of 
faces of Ci and C2. Furthermore, using the notation of Proposition 2.6, it is easy to see that 
dcixC2{^ij^2) = dci{xi) + dc2{x2). Hence, it follows from the characterization (15) that the 
intrinsic volumes of Ci x C2 arise as the convolution of the intrinsic volumes of Ci and C2. □ 
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Example 2.10. The calculation rule from Corollary 2.9 provides an easy way to compute 
the intrinsic volumes of the positive orthant: For n = 1 we have V{R+) = (i, |). As Rl = 
M_l_ X ... X M+ we may apply (18) to get 

V(W^) = (1, 1) * ... * (1, 1) = ((-)/2-, (^)/2«, . . . , , 
i.e., we have Vj-i{Rl) = (p /2" for Q < j < n. 

In the following proposition wc give a well-known formula for the spherical intrinsic volumes 
of smooth caps, which goes back to H. Weyl [57]. (Sec [57], [14] or [3, Ch. 4] for proofs.) 

Proposition 2.11. Let K G /C™(<S'"^~^) and 1 < j < n — 1. Then the intrinsic volumes of K 
are given by 

Vj-iiK) = / a„_,_i(p)dM, (19) 

where M := dK denotes the boundary of K, and a^ip) denotes the kth elementary symmetric 

function in the principal curvatures of M . □ 

Example 2.12. Let K = B{z,/3), < /3 < 7r/2, be a circular cap. Then from (19) and from 
Example 2.2 wc get for 1 < j < n — 1, with M := dK, 

Vj-i{K) = - ^3 / h~^).cot{/3r-'-UM 

f''~'^'^ .sm{/3y-'-cos{Pr-^-\ 



as vol(M) = sin(^)"'~^ • On-2- Furthermore, from (11) and recalling V-i{K) = Vn-i{K), we 
get 

Vn-i{K) = sin(p)"-2 dp , F_i(if) = ^ • cos(p)"-2 dp . 

^n-l Jo Jo 

We will finish this section by stating formulas for the intrinsic volumes of the cone of posi- 
tive semidefinite matrices. These formulas have been developed in [3], and a paper containing 
this derivation is in preparation [5]. We include these formulas here so that wc may appropri- 
ately formulate Conjecture 2.18 about the orders of magnitude of these intrinsic volumes. 

For z G we denote the Vandermonde determinant by A(z) := ni<i<j<A:(^i ~ ^j)- ^"^^ 
<r < k the Vandermonde determinant can be decomposed into 

r-(k-r) k-r 

Aiz) = A{x) ■ A{y) . ^ (-1)^-^ . pr4z) , Pr^z) := ® 2/"') -11^1' 

e=o i=i 

with X := {zi, . . . , Zr), y := {zr+i, . . . , z^), an denoting the ^th elementary symmetric function, 
and 

1 { Xj^ Xj^ X\ Xy ^ 

, . . . , , , . . . , , . . . , , . . . , 

y\ v\ y2 y2 Vk-r Vk-r 



X <i$ y 

Lastly, we define for < r < fe, < £ < r{k — r), and x, y as above. 



J{k,r,e) := J e-^ • |A(x)| • |A(y)| ■ Pr^z) dz , 



and we set J{k,r,£) := if the above inequalities on r and £ are not satisfied. 
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Theorem 2.13. The intrinsic volumes of the cone of positive semidefinite matrices are given 
by 

where t{r) := ^^'^^^^ = dim(Sym''). 

See the Figures 1,2 for graphical displays of the (logarithms of the) intrinsic volumes 
of Sym^ for k = 4,5, 6. The values of the intrinsic volumes were obtained by numerically 
approximating the integrals J{k, r,i) via a Monte-Carlo method. This explains the observable 
inaccuracies for A; = 6. 



2.3 Estimating intrinsic volumes 

In this section we will describe some quantities, which are related to binomial coefficients and 
which provide a method to estimate the intrinsic volumes of self-dual cones. Recall that in (9) 
we encountered the volume of the unit sphere, which we denote by On-i = vol^-i 5"^^. We 
denote the volume of the nth unit ball Bn C by 



LOji := VOlnBn = 



On- 



7r2 



n 



^^2n for n > 0, and := 1 • 



(20) 



It is convenient to use the analytic extension of the binomial coefficient. For x > — 1 and 
— 1 < y < X -|- 1 we denote := Y'{y+i)j\x-y+i) • particular, for half-integers, i.e., for 
n,m E 7., —1 < m < n + 1, we have 



n/2 
m/2 



Y^ in+2 _ n— m+2 ^ 



(21) 



Note that (, ) ^ ^n/2^ (("-2)/2) 

\(n-m)/2/ Vm/2/ V -1/2 / y^.r(lt±i) 

Besides the binomial coefficient we also use the flag coefficients [^] as defined in [34]. 
These are given by 



^•r(^) 



(22) 



2 J ^ y 2 I 

Note that = [^] . The following proposition provides some useful identities involving 

the binomial and the flag coefficients. 



Proposition 2.14. 1. For n,m eN, n > m, 



¥-r(^) •r(^^: 



r(^)-r( 2 

2. For n,m eN, n>m, 



n/2 
m/2 



(23) 



n/2 
m/2 



2 -On 



(24) 
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In particular, 



n 



OJn 



m ^n—m 



2-a 



n 



ITT' J ' ^n—m 



n/2 
m/2 



(25) 



3. For n, m — >■ oo such that also (n — m) — >■ oo, 



TT lm{n — m) f 

m/2 J ' 



2 V n 

where the symbol ~ means that the quotient of the two sides tends to one. 



(26) 



Proof. The first equation in (23) follows from applying the duplication formula of the T- 

The 



r(n+l) 



function r{2x) = ^ • 2^-'^ • r(x) ■ r(x + i) to the term (^) = r(n^+i).r(n-m+i) ■ 
remaining equations follow by plugging in the definitions of the corresponding quantities. As 
for the asymptotics stated in (26), we compute 



TT / m{n — m) f n/2 
m/2 



n 



m{n — m) 



m + 1) 



n 



r(f + i)-r(^^ + i) 



■•r(^). 



r(^^ 



r(^) •r( 



n—m+l ^ 



2 ' \l 2 ^ \ 2 ' 

where we have used the asymptotics ^/x ■ T{x) ^ V{x + 5) for x — >■ 00. 



□ 



With this notation we can revisit Example 2.12 and improve the formulas for the intrinsic 
volumes of circular caps. 

Example 2.15. Let K = B{z,P), < /3 < ir/2, be a circular cap. Then from Example 2.12 
we have for 1 < j < n — 1 



n-2 



Oj-i ■ On-j-1 



n 



■ sin(/3)^-i • cos(/3) 



(25) f{n - 2) /2\ sin(/3)J-i ■ cos(/3)"-J-i 



Furthermore, we have for j = n 

(11) O. " rl^ 



v r • ( \n-2, (21) An-2)/2\ n-1 ^ 



and similarly for j = 



n-2)/2\ n-1 fl'^ . . , „_2 ^ 

-1/2 )-^- 1 --(pr'dp. 
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For Lorentz caps, i.e., /3 = |, we have sin(|) = cos(|) = and the formulas for 
Vj-i{£,'"), which we denote by fj{n) for convenience, simphfy to 

/(n-2)/2x 

fjin) := Vj^.in = -^^^ , for 1 < j < n - 1 , (27) 
/o(n) := fn{n) := V.,{n = V^-iin = • 2("-2)/2 ■ {n - I) ■ £ sin(p)"-2 dp 



\ -1/2 } 
2"/2 



2Fi(l,i;^;-l) 



/(n-2)/2\ 

^ -1/2 / ^ 

t;^ — , tor n — >■ OO 

2n/2 ' 



where 2-^1 denotes the ordinary hypergeometric function (cf. [1, Ch. 15]). Note that the 
sequence /(n) is symmetric, i.e., /„_j(n) = fj{n). 

It is convenient to estimate the intrinsic volumes of self-dual cones by means of the intrinsic 
volumes of the Lorentz cones. We therefore define for a self-dual cone C C M" the excess over 
the Lorentz cone v{C) as 

v{C) := min . (28) 

0<i<n fj{n) 

In other words, v{C) is the smallest constant such that the inequality V{C) < v{C) ■ f{n) is 
satisfied componentwise. The constant v{C) will come in handy for the estimations, as we 
will see in Section 3. 

Proposition 2.16. For the positive orthant and for the Lorentz cone we have 



< V2 , !;(/:") = 1 . 



Proof. The claim about the Lorentz cone follows trivially from the definition of v and /. As 
for the positive orthant, we use Example 2.10 and compute for 1 < j < — 1 

v,.,{Ri) _ © v^.r(^).r(^).r(i±i).r(^) i 



/^■H ((J:?)/') • 2"/' r(i±i)-r(i±^).r(!^).r(^^)-r(f) 2^2 

(21) r- r( — ) [j/2) r- r( — ) l„/J /- 

where the last inequality is easily checked with a computer algebra system. As for j G {0, n}, 

we have 



< 1 



fn{n) foin) 2-/2 r(f).2Fi(l,i;^;-l) 

where again the last inequality is easily checked with a computer algebra system. □ 

An important construction in convex optimization is taking direct products. It is thus 
desirable to have an easy estimate for v{Ci x C2). Unfortunately, this is more challenging 
than it might appear at first sight. The following conjecture is supported both by geometric 
reasonings and computer-aided experiments. 
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2 4 6 8 10 5 10 15 5 10 15 20 

(a) A; = 4 (b) fe = 5 (c) fe = 6 



Figure 1: The intrinsic volumes of Sym^ (solid) and the upper bound 2 • f(t{k)) (dashed), 
illustrating Conjecture 2.18. 



2 4 6 8 10 5 10 15 5 10 15 20 




(a) fc = 4 (b) A; = 5 (c) fe = 6 



Figure 2: The logarithms of the intrinsic volumes of Sym^ (solid) and the logarithm of the 
upper bound 2 • f{t{k)) (dashed), illustrating Conjecture 2.18. 

Conjecture 2.17. For ni,n2 G N we have 

/(ni)*/(n2) <2-/(ni + n2) . (29) 

In particular, for closed convex cones Ci, C2 we have v{Ci x C2) < 2 • v{Ci) ■ f (C2), and 
thus 

v(>C"i X ... X £"'=)< 2^=-^ . (30) 

For the cone of positive semidefinite matrices we put up the following conjecture about 
the order of magnitude of its intrinsic volumes. 

Conjecture 2.18. For the cone Sym^ of positive semidefinite matrices we have 

u(Sym^) < 2 . (31) 
See Figure 1 and Figure 2 for an illustration of this conjecture. 
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The constant 2 in (29) and (31) is in both cases probably not optimal. In fact, concern- 
ing (29), computations suggest that f{ni)*f(n2) = (l + o(l)) • /(?^l +?t-2) for ni, 77-2 — )• 00. Our 
conjectures about the excess over the Lorentz cone are succinctly summarized in the table of 
Theorem 1.4. 

Concluding this section, we point out that the statistical behavior of intrinsic volumes is 

currently a field with a wealth of open questions. Note that Figure 1 shows an apparently uni- 
modal sequence, and Figure 2 suggests that it should actually be log-concave. We conjecture 
that this holds in general. 

Conjecture 2.19. For every K G 1C{S"'~^) the sequence of intrinsic volumes V{K) is log- 
concave, i.e., Vj{K)'^ > Vj-i{K) ■ Vj+i{K) for < J < n - 2. 

Remark 2.20. Using the formulas in Example 2.15 and the calculation rule from Corol- 
lary 2.9, one can show (cf. [3, Cor. 4.4.14]) that a direct product of circular cones has a 
log-concave sequence of intrinsic volumes. In particular, the intrinsic volumes of the positive 

orthant and of products of Lorentz cones are log-concavc. If Conjecture 2.19 is true, it could 
be considered as a spherical analog of the famous Alexandrov-Fenchel inequalities (cf. [52, 47]). 

3 A tube formula for the Grassmann manifold 

Before we state the main theorem we introduce the following notation for a rescaling of the 
volume of the tube of radius a around an i-dimensional subsphere (cf. (10)) 

where r := tan a. 

In the following, we fix a regular cone C C and set T>m ■= T>m{C), Vm ■= T^m{C), and 
Tjjn '■= ^m.{C), cf. Definition 1.1. Furthermore, we define the tube around as 

T{^m, a) := {W e GTn,m \ € : dg{W, W) < a} , (33) 

where dg denotes the geodesic distance in Grn,m (see [4] for more details). We also define the 
primal and the dual tube around S^, respectively, as 

r^(S„,a) :=r(S„,a)nP^ , r°(S„, a) := r(S„„ a) n . (34) 

In the following theorem we will give an estimate of the relative volume rvol of the primal 
tube around S^. By the relative volume we mean the unique probability measure on the 
homogeneous space Gr^^^, which is invariant under the action of the orthogonal group (see 
Section 5 for more details on this). The proof of Theorem 3.1 below is postponed to Section 6. 

Theorem 3.1. Let C C 6e a regular cone. Then, for 1 < m < n — 1 and < a < |, 

n-2 

(35) 

i=0 
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Table 1: The coefficient matrices Dn^j^ = (<if,^)ij=o,...,n-2 for n = 7, m = 1, . . . , 6. 



where the constants dfj^ are defined for i + j + m = 1 (mod 2), < + < m — 1, and 



Q < i+i _ <n-m-l, via 



/ m— 1 \ / n—m—L \ 

\ i—j I m— 1 ) ■ 1 i+j m — l ) 
2 2 2 2 



(36) 



and d^J^ := otherwise. 



Remark 3.2. • It can be sliown (see [3] for details) that (35) in fact holds with equality 

if T(C n 5""^, a) is convex and if d^^ is replaced by (—1)^"^^ • d^?". This reveals 
how close is our estimate (35) to being sharp. 



The coefficients d^p satisfy the symmetry relations 



fn,n—m mm 

^i,n-2-j ~ "'ij ' 



^nrn nun 

"'n—2—i,n—2—j '-''ij 



(37) 



See Table 1 for a display of the coefficient matrix Dn,m = {dij")i,j=o,...,n-2- 



The following corollary will be obtained from Theorem 3.1 by using duality properties as 
described in [4] to obtain an estimate of the dual tube of E^, and by changing the summation 
over i,j as depicted in Figure 3. 

Corollary 3.3. Let C C M" be a regular cone and let dfp be defined as in (36). Then, for 
1 < m < n — 1 and < a < |, and using the convetion (^) := if £ < or £ > k, we have 

n-2 p/ m+i-2fc+l i p/ n-Tn-i+2fc+l) \ 

rvoir(S^(C),a) <8 - ^ K (C) • ^ ' ■ ^ '- 

i,k=0 2 J 



m — 1\ f n — m — 1 
i — k 



r(^) 

In,n-2-i{oi) ■ (38) 
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Proof. The involution l: GTn,m Gr„^„_^, W i->- W"*" is an isometry and thus preserves the 
geodesic distance. Moreover, l maps Vrn{C) to Vn-rn{C), and therefore i(Em(C)) = S„_„j(C'), 
cf. [4] . Hence, i induces a bijection between {Yirn{C) , a) and iT,n-m{C) , a) . So we get 

rvol r(S^(C), a) rvol r^(S^(C), a) + rvol r°(S^(C), a) 

= rvol r^(S„(C), a) + rvol r^(S„_^(C'), a) , 

From Theorem 3.1 we get, using the duality property Vj{C) = Vn-2-j{C) (cf. (17)) and the 

-3 «J 

Am{n — m) /n/2' 



symmetry relation (i"^_2"!,- = d""^ (cf. (37)), 



rvol r(E^(C),a) < 



n 



m/2 



n-2 
j=0 



n-2 
L i . 



n-2 



^d^™ •/„,,(«). (39) 



1=0 



Using Proposition 2.14, we get 

(n/2\ rn-2] /'"/2\ 

4m(n - m) \m/2) ' l j \ (23) 4m(n - m) \m/2) 



n 



n 



/(n-2)/2N 
^ i/2 i 



(21) 4m(n — m) 



p(n±2) 



n 



p(m±2^ .p^n^m+2^ 



r(ag).r(^) 
r(?) 



r(^) r(^) 
r(f) r(^^) 



Using this in (39), changing the summation via i n — 2 — i and j n — 2 — j, and taking 
into account (37), we thus get 



(39) T(i^) rf^l 

rvol r(S„(C), a) < 8 • J] • • y,(C) • 

i,j=0 '-^2) H 2 ) 



n-2 



mm 



4,i(a) 



(37) g . g r(i±?) r(^^) 

i,j=0 
n-2 



r(f) r(^) 



n-2-j 



(C) 



n-2 



■ <^i^™ ■ In,n-2-i{a) 



ET/ (r^ ^ 2 ) 2 ) 
yn-2-jy^) ■ p/m\ ■ -rr n-m \ 
i,j=0 ^{ 2 ) 

•i+j+m=l 
(mod 2) 



m — i 

2 2 



n — m — 1 

i+j _ m-l 
2 2 



(40) 

In,n-2-i{(x) ■ 



Here we interpret 



if ^ < or £ > A;, i.e., the above summation over i,j in fact 



only runs over the rectangle determined by the inequalities < + ^^^^ < m — 1 and 
< — ^^^^^Y^ < n — m — 1 (cf. Figure 3). As the summation runs only over those for which 
i+j+m = 1 mod 2, we may replace the summation over j by a summation over A; = ^ + ^^^- 
The above inequalities then transform into < k < m — 1 and 0<i — k<n — m — 1. So we 
get from (40) 



n-2 



rvoir(S„(C),a)= 8 • ^ F„ 

—m—l—i+2k 



p^ m+i-2fc+l ^ p^ w-m-i+2fc+l) ^ 



i,k=0 



m — 1 
k 



r(f) 



n — m 



r(^) 

^n,n-2-i(a) • 



□ 
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Before giving the proofs of the main Theorems 1.3 and 1.4 by estimating the tube for- 
mula (38), we will provide a number of technical estimates in the following lemma. In the 
proof and in the rest of this section we will mark estimates, which are easily checked with a 
computer algebra system, with the symbol SI. 

Lemma 3.4. Let i, k, £,111,11 & N with n>2 and 1 < m < n — 1. 

1. We have 

r(=^)<,/!^7!!i±iy . (41) 



r(f ) - V 2 V 2 

2. For < k < m — 1 and 0<i — k<n — m — 1 we have 

i—2k 

m + i-2k \ ^ ^ . 
■■ — ^ ■ (42) 

3. For < a < |, t := sin(Q;)~^, and n>2>, we have 

E^'^fn — 2\i^ , . 3 , 3 

( . \ ■ ■ In,n-2-i{a) < -, ift>n2, (43) 



i=Q 

n-2 



\ ■ In,n-2-i{oi) < exp(^-J--, if t > m . 



(44) 



Proof. (1) We make a case distinction by the parity of £. Using r(a; + 1) = a; • r(a;), we get 
for odd £ 



n(?-) 



1 , I 

2 



r(f^) -LJ-V2 y-2V 2 / -V2V2 

Using additionally T[x + |) < ^/x ■ T(x), we get for even £ 



-1 , , , 1 . , < 

< 



r(^^^) _ /m + 1 \ ri^)^fm + £Y^ 



r(f) ni 2 r(f) 

(2) As for the second estimate, we distinguish the cases i >2k and i < 2k. From < < 
m — 1 and 0<i — A:<n — m — Iwe get 

1< m — k + i — k <n. — 1, 
1 < n — {m — k + i — k) < n — 1 . 

For i > 2A; we thus get 

i-2fc 

m + i — 2k \ 2 , . i-2k i 

< (n — 1) 2 < n2 , 



n — m — i + 2k 
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and for i < 2k 



i — 2k 2k — i 

m + i — 2k \ 2 /n — m + 2A; — A 2 , 2fc-» » 
r = ; <(n-l) 2 < na 



n — m — I 



(3) The /-functions have been estimated in [14, Lemma 2.2] in the following way. Let 
e := sin(a) = |. For i < n — 2 



In,n-2-i{oi) = I cos(p)"-^~^ • sin(p)* dp < 
Jo 

and for i = n — 2, assuming n > 3, 



.1+1 



-n-l 



2(n - 2) 



With these estimates we get 

n-2 



XI ( • ) • • In,n-2-i{a) 
n-2 



/ n - 2\ . e^+i 



n — 2\ i_ i 1.4 



j_o V ^ / ^ + 1 VV2(n-2) n-l 

] (^^(n-2\ , , 
< £ • I > • 712 . -I — — . £" - . n 2 

= £•1 (1 + •£)""' + 1.4 -e^-^.n"^' 



• e" • ra 2 



_3 

For £ < n 2 we thus get 



2 . nt • /n,n-2-^(a) < £ . (1 + iV ' +1.4 • ni-") < 3 • £ 



i=0 



v 

<exp(l) 

Similarly we derive 

n-2 



E ') • ^n,n-2-.(«) < e . ((1 + £)-^ + ^ • £-^) 

For £ < ^ we thus get 

> . • /„,„_2-i(«) < £ • 1 + - + ^ ^ < £ • exp 



i=0 



□ 



20 



Proof of Theorem 1.3. If C is a regular cone then it is not a linear subspace. The intrinsic 
volumes arc therefore bounded by Vj{C) < ^ (cf. Remark 2.8). Using Lemma 3.4 we get from 
Corollary 3.3 

(38) p/ m+i-2fc+l \ J^/ n-m-^+2fc+l) ^ 

rvol r(S^, a) < 8 • ^ —m—i+2k—l 



r(f) r(^) 

• In,n-2-i{oi) 



.,k=0 

^m — 1\ fn — m — 1 
k )'\ i-k 

(41) ^- ^ fm + i-2k\~ fn-m-i + 2k\ ^ 

< 2 V"i(n-m)- { ^ •( ^ 

m — 1\ f n, — m — 1 



k I V i — A; 



• In,n-2-i{a) 

(42) ^ — ^L-^ i /m — l\ f n — m — I 

o — n z — n \ 



i=0 A;=0 

n— 2 /m— 1\ /n— m— 1\ /n— 2\ 



/ V i — k 



By Vandermonde's identity we have Y.l=o ("fc ) ' {'\-k ) = ("7 )• So we get 

■n-2 / _ 2\ 

rvol T(i;m, a) < 2 • V'"^(n - m) • ^ ( . ) • nz • /„^„_2_i(a) 



(43) / 1 3 

< 6 • V m(n — m) • - , if t > n 2 and n > 3 . 



As for the expectation of the logarithm of the Grassmann condition, we compute 

POD 

E[ln(^G(^))] = / Prob[ln(^G(^)) > s]ds 
Jo 

roo 

< 1.5 • ln(n) + r + / 6 • -\/m(n--ni) • exp(— s) 

^ln(n'V2)+r 

Jmiri — m) 
= 1.5 • ln(n) + r + 6 • — ^-t;^ • exp(— r) 

' 

<|a/6 , if n > 3 

< 1.5 • ln(n) + 1.5 , 

if we choose r := ^ In (|) . □ 

Before we finish this section with the proof of Theorem 1.4 we need yet another technical 
lemma. A sequence (a„) of nonnegative real numbers is called log-concave iff > fln-i ■ ^n+i 
for all n. See [52] for a survey on log-concave sequences and their appearances in diverse areas 
of mathematics. 

Lemma 3.5. For n > 2 and 1 < m < n - 1 let gm{n) := ^^y^^^^y?^- 

1. The sequence (5m(^))n is log-concave, i.e., gminY > 5m(^~l)"9'm(^ + l) forn > m + 2. 
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2. For fixed m> 8 we have max{^^(n) | n > m} = max{^^(2m + k) \ k E {5, 6, 7}}. 

3. For m> 8 we have 



gm{n) <2.5-y/m. (45) 

Proof. (1) We have 

ff^(n)^ _ r(§)^ r(^^) TC-^) 

gmin-l)-gmin + l) r(^)2 r(^^) r(^) " 

In order to show that this expression is > 1 we use induction on m. For m = this is triviahy 
true, and for m = 1 this is easily checked with a computer algebra system. For m > 2 we 
have, using r(x + 1) = .x ■ T{x), 

r(§)2 r(^^^^) r(^^±^) 



r(rwn)2 r(n^) r(n±l) 



2 

(n-m)2 r(f)2 r(n-l-|rn-2)^ „+i-(„,-2) . 



(n - 1 - m) • (n + 1 -m) p( n-(m-2) ^2 r(^) r(^) 



> 1 . 



>i 



>1 by ind. hyp. 



(2) As the sequence {gm{n))n is log-concave and positive, it follows that it is unimodal. 
This means that there exists an index N such that gm{n — 1) < gmin) for all n < N, and 
gm(,n) > gm{n + 1) for all n > iV (cf. [52]). Moreover, for m > 8 we have A'' G {2m + k \ k & 
{5,6,7}}, as 

gm{2m + 4) _ Ti'-^) r(^) ^2 J ^ 



gm{2m + 5) r{^) r(^^) exp(i^) 
5^(2m + 7) _r(2zr^) r(^rf8) ^ 



gm{2m + 8) Ti^) r(^) exp(^) 
(3) For fixed k, the following asymptotics is easily verified: 

gm{2m + k) m-^oo 4-exp(2) 



> 1 , for m > 8 . 



m V27r 



< 12 



In particular, it follows by (2) that for m > 8 we have an asymptotic estimate of gmin) = 
0{y/m). More precisely, it is straightforward to check that for k G {5,6,7} and m > 8 we 
have gm{2m + k) < 2.5 • s/m. It follows by (2) that for m > 8 we have gm{n) < 2.5 • ^/m. □ 



Proof of Theorem 1.4- We will estimate the intrinsic volumes of C via V^_i(C) < v{C) ■ fj{n), 
where fj{n) = Vj_i(>C") and v{C) denotes the excess over the Lorentz cone introduced in 
Section 2.3. Note that for j = m + i — 2k with 1 < j < n — 1, we get from (27) using (21) 

f ^ y{n-m-i+2k-l)/2J ^ U2J , . 

Jn-{m+i-2k){"') ^71/2 p^ n-m-i+2fc+l p|- m+i-2fc+l ■) ' ^ 
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We thus obtain from Corollary 3.3 

(38) p/ m+t-2fc+l \ p/ n-r»--i+2fc+l) \ 

rvol r(S„,a) < 8 • 5^ Vn-m-i+2k-i{C) ■ I ■ ' 

i,k=0 n 2 ) 

(m — l\ (n — m — l\ 



(46) r(~) 

^ 8 • V{C) ■ ^;7^^-^;7^^^^;w^ • In,n-2-i{a) • ^ T 

r(-) "^^ fn — 2\ 

8 • 'f^(C) • r(m^ r(n-my^^ " In,n-2~i{a) " ^ J 



A; / V i — k 



Using the notation from Lemma 3.5, this implies 

1 (45) ^ 1 

rvol Ti^m, a) < 8 • v{C) ■ gm{n) ■ - < 20 • v{C) • • - , for t > m > 8 . 

Analogous to the proof of Theorem 1.3 we estimate the expectation of the logarithm of 
the Grassmann condition. Defining v{C) := max{t;(C), 1}, so that in particular \nv{C) > 0, 
we get 

E [ln(^G(A))] = / Prob[ln(^G(^)) > s] ds 
Jo 

POO 

< ln(m) + ln(u(C)) + r + / 20 • v{C) ■ yfm ■ exp(-s) ds 

Jlnfml+lnfCl'CD+r 

= ln(m) + ln({;(C)) + r 



'ln(m)+ln({;(C))+r 

20-?;(C) 



fm ■ v{C) ■ exp(r) 
< ln(m) + ln({)(C)) + 3 , if m > 8 , 

where the last inequality follows from choosing r := In □ 

4 Twisted characteristic polynomials 

This section is about a certain subspace-dependent version of the characteristic polynomial of 
a linear operator. This polynomial will appear in the computation of the volume of the tube 
around the set in the Grassmann manifold. Nevertheless, we will formulate what we call 
the twisted characteristic polynomial in a general context. The proof of the main result of this 
section is elementary and only involves basic linear algebra. 

Let (p be an endomorphism on a A;- dimensional euclidean vector space V . We denote by 
crj{ip), < j < k, the coefficients of the characteristic polynomial of ip (up to sign). More 
precisely, 



det{ip - t ■ idy) = Xl(-l)''"' • ^i(^) • 

i=0 
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Note that we have (Tfc((^) = det{ip), cro((/?) = 1, and cri{(f) = trace(<^). In the following we 
denote by Gr(F, £) the set of all ^-dimensional linear subspaces of V. 

Definition 4.1. Let Y G Gr(y,l), and let Ily and Ilyi denote the orthogonal projections 
onto Y and Y^, respectively. The twisted characteristic polynomial chy and the positive 
twisted characteristic polynomial chy of ip with respect to Y are defined via 

chy(v?,t) := dot {if - {t-Uy-j ■ Hyx)) • t'''^ , 
ch+{<f, t) := det + (t • Hy + I • Hyx)) • t'"'^ . 

Furthermore, we denote by (py : Y ^ Y the restriction of (p, i.e., <PY{y) ■= Ily ((^(y)), and we 
use the notation dety(99) := dct((/3y). 

Note that for ^ = A; we get chyi^, t) = det((/3 — t • idy), the usual characteristic polynomial, 
whereas for £ = we get cho((^, t) = det(t • ip + idy). We will see below that 

chy (^, 0) = ch+ (^, 0) = dety {^) . (47) 

Theorem 4.2. Let V he a k-dimensional euclidean vector space and let (p be an endomorphism 
on V. IfYG Gi{V,£) is chosen uniformly at random, then 

E[dety{^)] = ^ ■ aeiip) . (48) 
Moreover, the expectation of the (positive) twisted characteristic polynomial is given by 

k k 

E [chy (v., t)\ = J2 dij ■ ('k-jiv) ■ i'-* , I [ch+(^, t)] = J2 \dij\ ■ ^^k-M ■ t^'' , (49) 

i,j=0 i,j=0 

where the coefficients dij are given for i + j + i = mod 2 and < ^ + | < -^^ < ^ ~ i ^ 
k — i, by 

and dij := otherwise. (Note that \dij\ = d\^'^'^'^^ in the notation of (36) .j 

If ip is positive semidefinite, then for t > we have |ch-y(y3, t)| < chy((^,f) and 



k 

chy{ip,t)\ < XI I I ■ '^'^-J ('^^ ■ 

i,j=0 



xk 



We will now express the twisted characteristic polynomial in coordinates. For A G 
and < £ < k we use the notation 

ch^(^,t) := chT^i^o{(p,t) , d4{A,t) := ch+ ^^((^,4) , 

where : M*^ — >■ M'^, x i->- Ax. Note that if A has the block decomposition A= ^ ^ , where 
Ai G M^^^, and the other blocks accordingly, then 

Ai-tii A2 \ ,fAi + tIi A2 ^ 



ch,(A«) = det( - ^^^^',] , ch;(A,<) = det - ^^^^',] . (51) 
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From this description we get for F = X the identity chy(<^,0) = ch+((^,0) = det(Ai) = 
dety ((/?). This proves (47). 

If lb £ Gik,i is fixed and Q € 0{k) is chosen uniformly at random, then the induced 
probabiUty distribution on Gr^^^ via Q QYq is the uniform distribution. Prom this it is 
straightforward to check that if Q G 0{k) and Y G Gr^^^ are chosen uniformly at random, 
then we have 

E[chYi<p,t)]=E[che{Q^AQ,t)], E [ch+ {ip, t)] =E[ch+ {Q'^ A Q,t)] . (52) 

These equalities allow to prove Theorem 4.2 with basic matrix calculus. The main idea is to 
use the multilinearity of the determinant and the invariance of the coefficients aj{A) under 
similarity transformations to show that the coefficients of the expectation of the (positive) 
twisted characteristic polynomial arc linear combinations of the a.j(A). Wc may then compute 
the exact linear combinations by choosing particularly nice matrices A, namely scalar multiples 
of the identity matrix. 

In the following we use the notation [k] := {1, .... k}, and wc denote by ('^') the set of all 

i-element subsets of [k]. For A G M*^^^' and J G (^f), we denote by pmj(A) the Jth principal 
minor of A, i.e., pmj(74) := det(Aj), where Aj denotes the submatrix of A which results 
from selecting the rows and columns of A whose indices lie in J. It is well-known (cf. for 
example [32, Thm. 1.2.12]) that ai{A) is the sum of all principal minors of A of size i, i.e., 

aiiA)= prajiA). (53) 

Mm 

For < i < k the ^th leading principal minor is the principal minor for J = [£]. We denote it 
by Ipm^(^) := pm[^](A). Note that 

Ipm^(A) = det^e^oi^p) . 

Lemma 4.3. Let A G M*'^'^, and let Q G 0{k) be chosen uniformly at random. Then, for 
J G ('^') , we have 

E [pm^(Q^ -A-Q)] =E [lpm^(Q^ ■ A ■ Q)] = -J- • ae{A) . (54) 

Proof. For the first equality let J = {j'l, . . . ,j£}, ji < . . . < ji, and let tt be any permutation 
of [k] such that Tr{i) = ji for all i = I, . . . , L If Af^ denotes the permutation matrix according 
to TT, i.e., Mt^-Ci = e^(j), then Aj = {MJ -A-M^)^;^^, and therefore pmj(^) = lpm^(Mj-yl-M;r). 
This implies 

E [pm^(Q^ ■A-Q)]=E [lpm^(Mj ■ Q'^ ■ A ■ Q ■ M^)] = E [lpm^(Q^ -A-Q)] , 

where we have used the fact that right multiplication by the fixed element leaves the 
uniform distribution on 0{k) invariant. This also implies 

E[lpm,{Q^.A.Q)]=^. ^ E[pmj{Q^.A.Q)] = ^.E\ ^ pmAQ^-A-Q) 



25 



Proof of Theorem 4-2. The statement (48) follows from (52) and (54). For the equalities 
in (49) it remains to show that for Q e 0{k) uniformly at random 

k 

E [cheiQ^'AQ, t)] = dij ■ (Tk-M) ■ t''~' > (55) 
^ i,j=o 

k 

E [ch+{Q^AQ,t)] = J2 ■ <^k-j{A) ■ ■ (56) 

Q ■ ■ n 

1,3=0 

By multilinearity we may write the determinant of a matrix with columns vi, . . . , Vi-i,Vi + 
a ■ Si, Vi+i, . . . , in the form 

det(i;i, . . . ,Vi-i,Vi + a-ei,Vi+i, ...,Vn) = det{vi, . . . , v„) + a ■ det('yi, . . . , ■Uj_i, e^, Vj+i, . . .,Vn)- 
Using this repeatedly, wc generally obtain 

det(A + diag(ai,...,an)) = ^ pmj(A) • Oj . 

jc[k] ie[fe]\J 

As we have (cf. (51)) 

ch^(^, t) = t''-^ ■ det{A + diag(-t, . . . , -i, . . . , l/t)) , 
ch+(^,t) = 1^"^ ■diet{A + d\&g{t,...,t,\/t,...,l/t)) , 

we may expand the (positive) twisted characteristic polynomial to obtain 

ch,(A, t)=Y, i-lf^^'^ ■ v^j{A) . f-^') , ch+iA, t)=J2 Pmj(^) • t"'^'^ ' (57) 

JQ[k] JQ[k] 

where Ci,C2: 2''^] — >■ N are some integer valued functions on the power set of [A;]. Averaging 
the twisted characteristic polynomial thus yields 

E [ch,(Q^ .A-Q,t)]=J2 i-^y^'^ ■ ^ [pmj(Q^ -A-Q)]- e-^') 

JC[fe] 

JC[fe] \\J\) i,j=0 

for some rational constants dij. To compute these constants, let us consider the matrices 
A = s ■ Ik- For this choice of A we have ch.£{s ■ 1^, t) = (s — tY ■ {1 + s ■ t)^^^ . As ak-j{s ■ 1^) = 
Q ■ s^-^, and Q'^ ■ s- Ik- Q = s- Ik, we get 

{s-t)'-{l + s- tf-^ = cMs ■ Ik, t)=E [cMQ^'^ ■s-h-Q,t)] = Y, 4 • ■ ■ i^'' ■ 

^ i,j=0 

We expand the first term in order to make a comparison of the coefHcients to get the dij . We 
have 
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Figure 3: Illustration of the change of summation in (58) {k = 8, £ = 5). 




(mod 2) 



where again we interpret (^) = if m < or rn > n, i.e., the above summation over i,j 
in fact only runs over the rectangle determined by the inequalities < + ^ < ^ and 
< A: - ^±|±^ <k-i. See Fi gure 3 for an illustration of the change of summation. Note that 
the reverse substitution is given by A = ^—^^ and ji = k — 

Comparing the coefficients of the two expressions of (s — f)^ • (1 + s • t)^"^ reveals that 
indeed dij = dij as defined in (50). This shows the equality in (55). 

The equality in (56) is shown analogously with the observation 

c\4{s-Ik,t) = {s + tf-{l + s-tf-^. 

As for the last claim in Theorem 4.2, note that for positive semidefinite A every principal 
minor is nonnegative, i.e., pmj(^) > for all J C [k\. Therefore, if t > 0, we get from (57) 

|ch,(^,t)| I ^ (-l)^^(^) •pm^(A) < ^ pm^(^) -i^^^-^) ch+(A,t) . 

JC[fe] JC[fe] 

This finishes the proof of Theorem 4.2. □ 

5 Preliminaries from Riemannian geometry 

We refer to [19], [8], or [15, Ch. 1], for background on Riemannian geometry. In the following 
subsections we will focus on three different subjects. First, we will state the smooth coarea 
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formula, which is the basis for all subsequent volume computations. Second, we will describe 
the Grassmann manifold as a quotient of the orthogonal group. This viewpoint is a great 
help in doing explicit calculations in the Grassmann manifold. The third and final topic is 
a description of the (orthonormal) frame bundle and the Grassmann bundle of a manifold. 
These bundles form the geometric basis for the computation of the tube formula in Section 6. 



5.1 Coetrea formula 

An important tool in our computations will be the smooth coarea formula. Before we can 
state this, we need to define the normal determinant of a linear operator. 

If A: V ^ W is a linear operator between euclidean vector spaces V and W, then we 
define the normal determinant of A as 

ndet(yl) := | det(yl|ker(A)-L)| , 

where ^|kcr(A)^ • ker(A)-'- — t- im(^) denotes the restriction of A to the orthogonal complement 
of the kernel of A. Obviously, if ^ is a bijective linear operator then ndet(^) = | det(74)|, so the 
normal determinant provides a natural generalization of the absolute value of the determinant. 

Theorem 5.1. Let (p: Aii — )■ A^2 be a smooth surjective map between Riemannian mani- 
folds M.i,M.2- Then for any f : Mi ^ M that is integrable w.r.t. dMi we have 

If additionally dim Ali = dim.M.2, then 

Yo\M2:= J dM2< J #(p~^{q)dM2= J ndet{Dp(p) dMi , (60) 

M2 qEM2 pEMi 

where ij^f~^{q) denotes the number of elements in the fiber ip~^{q). □ 

The inner integral in (59) over the fiber ip^^[q) is well-defined for almost all q G Ai2- This 
follows from Sard's lemma (cf. for example [50, Thm. 3-14]), which implies that almost all 
q € M2 are regular values, i.e., the differential Dpip has full rank for all p G ip~^{q). The 
fibers ip~^{q) of regular values q are smooth submanifolds of A^i and therefore the integral 
over ip~^{q) is well-defined. One calls iidet{Dpip) the Normal Jacobian of 9? at p. 

See [37, 3.8] or [22, 3.2.11] for proofs of the coarea formula with M.i,M.2 being subman- 
ifolds of euclidean space. See [33, Appendix] for a proof of the coarea formula in the above 
stated form. 



5.2 Orthogonal group and Grassmann manifold 

The Lie algebra of the orthogonal group is given by 

Ti^O{n) = Skew„ := {U G M"""" | = -U} . 
The tangent space of an element Q G 0{n) is thus given by TQO{n) = Q ■ Skew, 



28 



As for the Riemannian metric on 0{n), it is convenient to scale the (euchdean) metric 
induced from W^^"' by a factor of ^. The Riemannian metric (., .)q on TQO{n) = Q ■ Skew„ 
is thus given by 

{QUuQU2)Q=^^-tT{Ul ■U2) . (61) 
for Ui, U2 G Skewn- Observe that we have a canonical basis for Skewn given by {Eij — E, 



1 < J < ^ < n}-, where Eij denotes the {i,j)th elementary matrix, i.e., the matrix whose 
entries are zero everywhere except for the (i,j)th entry, which is 1. This basis is orthogonal 
and by the choice of the scaling factor it is also orthonormal. 

The exponential map of a compact Riemannian manifold is a map exp : TM. — >■ M. , 
where TM. denotes the tangent bundle of M., i.e., the disjoint union of all tangent spaces. For 
p & M and v G TpM., the map ^{t) := expp(tf ) is the constant speed geodesic with 7(0) = p 
in direction 7(0) = v. We have dg{'~f{to),^{ti)) = |to — ti\ • \\v\\ for \to — ti\ sufficiently small. 
Moreover, we have ^^(7(^0), 7(^1)) < \to — ti\ ■ \\v\\ for all to,ti. 

For the orthogonal group 0(n), the exponential map expg : TQO{n) —)• 0(n) at Q G 0{n) 
is given by the usual matrix exponential, i.e., for U G Skew„ we have ex.pQ{QU) = Q ■ = 
Q ■ Yl'k=o 1^- example, for 





( 








-1 \ 


N := 













V 


1 





/ 



(62) 



the exponential map at Q in direction QN is given by 





( 


COS p 





— sinp 


expQ(/3 • QN) = Q-Qp, Qp:= 





In-2 







V 


sinp 





COS p 



(63) 



Note that ^Qp = Qp-N. 

It can be shown (cf. for example [8, §VII.6]) that for sufficiently small r > the restriction 
of the exponential map at Q to an open ball BQ{r) of radius r around the origin in TQO{n), 
is a diffeomorphism between Bqir) and its image. This statement in fact holds for arbitrary 
Riemannian manifolds Ai. Furthermore, for arbitrary r > we have (cf. [8, §VII.6]) for p G M 
and Bp(r) C TpM. the closed ball of radius r around the origin in TpM.^ 

expp(5p(r)) = {qeM \ dg{p, g) < r} , 

where dg denotes the geodesic distance in M.. Moreover, if H C is a compact hypersurface 
in M. with unit normal vector field u, then we can characterize the tube of radius r around % 
as follows: 



r(H,r) := |g G I Tmadg{p,q) < r} = |J expp([-r,r] • v{p)) . 



(64) 



(For a proof see for example [51, Addendum to Chap. 9, proof of Thm. 20].) 
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The Grassmann manifold Gin^m, 1 < m < n, is defined as the set of all m-dimensional 
subspaces of R". We can describe Gr^^^ as the quotient 0{n)/H := {QH | Q G 0{n)} of 
0{n) by the subgroup 

ii" := I (^^ Qf G 0(m) , Q" G 0(n - m)| ~ 0(m) X 0(n - m) 

in the following way: 



0{n) 



0{n)/H 



n 



Gr„ 5 



, with n: g ^ Q(M™ X 0) . 



(65) 



In other words, the identification of Gr„^m with the quotient 0{n)/H amounts to identifying 
an m-dimensional linear subspace W of with the set of all orthogonal matrices whose first 
m columns span W. Note that this identification endows Gr^^^ with a topology, namely the 
quotient topology on 0{n)/H. 

In the following paragraphs we will give a concrete description of the tangent space 
n m for ^ £ Gr„ „j, and we will describe the Riemannian metric and the exponen- 
tial map on Gr„^m- Note that the natural action of 0{n) on M" induces a corresponding 
action on Gr^^^. It can be shown that, up to scaling, there exists a unique Riemannian metric 
on Gvn^rn-, which is 0(n)-invariant. This is the Riemannian metric, that we will describe next 
explicitly. 

Let G G Tn,rn sjvA let Q G 11 ^(VF), i.c, the first m columns of Q span W. The coset 
QH is a submanifold of 0(n), and its tangent space at Q is given by 



TqQH = Q ■ Ti^H = Q ■ 



U' 
U" 



U' G Skew^ , U" G Skew„_, 



Note that the orthogonal complement (TqQH)-^ of TqQH in TQO{n) is given by 



{TqQH)^ = Q ■ Skewn , where Skew„ := 



-X^ 
X 



(66) 



It can be shown (cf. [30, Lemma II.4.1]) that there exists an open ball B around the origin 
in TgOin) such that the restriction of 11 o expQ to i? H {TqQH)^ is a diffcomorphism onto 
an open neighborhood olW = Tl{Q) in Gr^^m. The derivative of 11 o expQ therefore yields a 
linear isomorphism 



{TqQH)- 



Tw Gr„ i 



QU ^i:=D{Jloe^VQ){QU) , 



where U G Skcw„. 

The pair (Q, U) with U G Skcw„ thus defines the tangent vector ^ := Dili o expQ((5^7)) 
in Ty/ Gin,m- However, we may also represent the subspace II/' by a group element Qh, for any 
h e H. A small computation shows that the pair {Qh,h~^Uh) represents the same tangent 
vector ^. More generally, if we define an equivalence relation ~ on 0{n) x H via 



{Q,U)^{Q',U') 



3heH: {Q',U') = {Qh,h-^Uh) , 
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then it follows that the tangent vector ^ G T\y GTn,m can be identified with the equivalence 
class [Q, U] := {{Qh, h~^Uh) \ h e H} e (0(n) xH)/ ^. We have thus obtained the following 
model for Grn,m: 

Tw GTn,m ^{[Q,U]\Q& ir\W) , U e Sk^} . 

Note that the derivative DqH of 11 at Q is given by 

DQn{QU) = [Q,TriU)], (67) 

where tt: Skew — t- Skew„ denotes the orthogonal projection. 

We can now define the Riemannian metric on Gr^^^ by setting 

{[Q,U,], [Q,U2]) := {Ui,U2)iJ=^ l-tviUf ■U.2) , (68) 

for Ui,U2 G Skew„ (this is clearly a well-defined 0(n)-invariant Riemannian metric). This 
way, 11: 0(n) — >■ Gr^^^, becomes a Riemannian submersion, i.e., for every Q G 0(n) the 
restriction of D expg to the orthogonal complement of its kernel is an isometry. Using the 
above description of DH and elementary properties of geodesies (cf. for example [19, Sec. 3.2]), 
one can show that the exponential map on Gr^^^ is given by 

exp^(0 = n(expQ(QC/)) = n(Q • e^) , ^ = [Q, U] G Tw Gr„,^ . (69) 

In Definition 1.1 we assigned to a regular cone C C and 1 < m < n — 1 the compact 
subsets Vm = 'Pm(C) and Vm = X'm(C') of the Grassmann manifold GTn,m- In Section 6 we 
will prove that S„j = Vm Pi T>m is a compact hypersurface in Gr„^m- Furthermore, wc will 
identify a unit normal vector field jy in which points into the component i.e., for 
a > small enough we have expp^(a • u) G Vm for all W G S^. In this situation, we can 
refine the characterization (64) of the tube around to a characterization of the primal tube 
T^(S^, a) = T(S^, a) n Vm as follows: 

r^(E^, a) = ( U expvK([0, a] ■ u{W))) n Vm ■ (70) 

The inclusion 'D' is clear by (64). As for the other inclusion, let W' G Vm with a := 
dg{W' ,T:m) > 0. Prom (64) we get that a minimum length geodesic between some W G 
and W' has the form 7(p) = exp^iS ■ p ■ i'{W)) with < p < a, and S G {1,-1}. It 
is sufficient to show that 6 = 1. By way of contradiction, assume that 6 = —1. Then 
W' = 7(a) = expp^(— a • iy{W)) G Vm- Since 7(p) G T>m for sufficiently small < p < a, 
there exists p < ao < a such that 7(0:0) € 'Pm H = ^m- Therefore, ^^(W^',!;^) < 
dp(7(a), 7(ao)) < a — ao < which is a contradiction. 

We may use the smooth coarea formula (59) to compute the volumes of the orthogonal 
group and of the Grassmann manifold. Note that the map ip: 0{n) 5""^, Q ^ Qei, where 
ei G M" denotes the first canonical basis vector, is a Riemannian submersion. In particular, we 
have ndet(Z?Q(/?) = 1 for all Q G 0(n). Moreover, each fiber of (p is isometric to 0{n — 1). So 
an application of the smooth coarea formula (59) shows that vol 0(n) = On-i " vol 0(n — 1). 
By induction on n we thus get volO(ra) = nr=o^^«- volume of GiCn,m, recall 

that by definition of the Riemannian metric on Gr„_^, we have a Riemannian submersion 
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n: 0(n) ^ Gr 

n^m whose fibers are isometric to the direct product 0(rn) x 0{n — m). The 
smooth coarea formula (59) thus imphes 

volGr vol 0(n) ^ Il7=n-m<^i 

"'"^ vol 0(m) • vol 0(n - m) J^^-i d 

For later use, we also compute 

V0lGr„_2,m-l _ Ili=n-m-l Yli=0^ _ ^m-l ' On-m-l (20) m(n - m) 



volGr„,^ l\T=fO^ \Yl=n-^0, On-2-On-i u 0„_2 • o;, 

(24) m(n — m) fn/2\ 1 



n \m/2) On-2 

5.3 Frame bundle and Grassmann bundle 



(71) 



In the remainder of this section we will discuss some elementary properties of frame bundles 
and Grassmann bundles. Note that the only difference between R"" and a general n-dimensional 
euclidean space is that the space M" has a distinguished orthonormal basis. 

Let F be a /c-dimensional euclidean vector space. We denote the set of orthonormal bases 
(frames) in V by F{V), i.e., 

F{V) = {{qi, ...,qk)eV^\ Vi, j : {q^, qj) = 5ij} . 

Note that one can naturally identify F{W^) with the orthogonal group 0{n). In particular, 
F(y) has the structure of a smooth (Riemannian) manifold. We recall the ^th Grassmann 
manifold Gv{V,£) = {Y C V \ Y an ^-dimensional linear subspace}, defined for 1 < £ < 
k. Again, since Gr(F, ^) 2± Grjt^^, the set Gr(F, ^) has a well-defined structure of a smooth 
(Riemannian) manifold. Note that we have the smooth surjective map F(y) — >■ Gr(y,£), 
{qi, . . . , qk) ^ lin{gi, . . . , qi}, whose definition is related to (65). 

Let M be a fe-dimensional smooth Riemannian manifold. The orthonormal frame bundle 
and the ^th Grassmann bundle, 1 < £ < k, over the manifold M are defined as 

F{M) = U M X FiTpM) , Gi{M,£) = \J {p} x Gi{TpM,£) . 

Using charts of M to connect the manifolds F{TpM), it is possible to show that F(M) has the 
structure of a smooth manifold. Moreover, F{M) ^ M is a fiber bundle (cf. for example [35, 
§1.5]). The same holds true for the Grassmann bundle Gt{M,£) — >■ M. Furthermore, the maps 
of the fibers F{TpM) Gv{TpM, £) may be combined to a smooth surjective (bundle) map 

n^:F{M)^Gv{M,£), Il^(j>,Q) = (pM{q^, . . . ^q^}) , (72) 

where =_(9i, ••• 

For [p, Q) G F{M) one has a natural decomposition of the tangent space of F{M) at {p, Q) 
into a vertical space and a horizontal space, 

TM)nM) = T^,Q)F{M) © r^,Q)F(M) . (73) 
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The vertical space is defined by 



%^Q)F{M) := T(p,g)(M X F{T,M)) . (74) 

For defining tlie horizontal space we consider the following construction. Let Q G XpM and 
let c: M — )■ M be such that c(0) = p and c(0) = C,. Furthermore, let Qt denote the parallel 
transport (cf. for example [8, Thm. VII.3.12]) of the frame Q along c at time t. Then the map 
c^: M — ^ F{M), t i-> {c{t),Qt), defines a smooth curve in F{M) and thus a tangent vector 
cir(O) G T^pQ-^F{M). It can be shown (cf. for example [35, §11.3]) that this tangent vector docs 
not depend on the specific choice of the curve c. Moreover, one can show that the induced 
map TpM T^p Q^F{M) sending each tangent vector ( to the above defined tangent vector 
of F{M) at (p, Q), is an injective linear map. The image of this injective linear map is defined 
as the horizontal space. One can then show the decomposition (73) of the tangent space of 
F{M) at (p, Q) into the direct sum of the vertical and the horizontal space. 

Analogous statements about the decomposition of the tangent spaces also hold for the 
Grassmann bundle Gt{M,£). Moreover, the derivative of the bundle map I!;,: F{M) 
Gr{M,i) certainly maps the vertical spaces of F{M) to the vertical spaces of Gi{V,C^. As for 
the horizontal spaces, note that the parallel transport of a linear subspace may be achieved by 
parallel transporting a basis for this subspace. This implies that also the horizontal spaces of 
F{M) are mapped to the horizontal spaces of Gtl{V, I) via the derivative of the bundle map n^. 

6 Establishing the tube formula 

The following theorem is the main result of this section. Before approaching its proof, we will 
use it to give a proof of Theorem 3.1. 

Theorem 6.1. Let K G /C="^(S'"-i), C := conc(A'), and let M := dK. 

1. The set Tim '■= Tm{C) is an oriented hypersurface ofGin,m- The map 

Gr(M,m-l)^S^, {p,Y)^Rp + Y (75) 
is a diffeomorphism with inverse given by $~^(VF) = {p,Y), where W D K = {p} and 

Y =p^nw. 

2. For W G T,m let v-^^W) G T\y GiCn,m denote the normal direction pointing inside V^iK), 
and let the maps Hm and * be defined via 

Hm : Sm ^ M ^: S„ X M ^ Gr„,^ (76) 

W ^ p , where W = {p} , (l^", t) ^ exp^;(/(arctan(-t) • v^^iW)) . 

Then the Normal Jacobians of Hm and are given by 

ndet{DwU^) = deiY{Wp)-\ | det(i^(t^,,)^')| = (1 + i2)-/2 . ^ (77) 

where {p,Y) = ^^(W), and Wp denotes the Weingarten map of M at p. 
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Proof of Theorem 3.1. Recall that the intrinsic volumes Vj{K) are continuous in K (cf. Propo- 
sition 2.4). Moreover, the set of smooth caps /C^™(S'"~^) lies dense in /C'^(S'"~"'^); in particular, 
it lies dense in the set of regular caps K7{S'^~'^) (cf. Section 2.1). Furthermore, it is straight- 
forward to check [3, Lemma 6.1.7] that the map 

/C(S"-^)^M, rvoir°(Sm(i^),a) 

is continuous and thus, by compactness of /C(S'"~^), uniformly continuous. Therefore, in order 
to show the inequality (35), we may assume w.l.o.g. that K G /C™(S'"^"'^). 

Let the notation be as in Theorem 6.1. As the unit normal field v^. of points into the 
component Vm{K), i-c, the direction —u^ points into the component Vra{K), we have for 
T := tan a 

*(E^ X [0, r]) D ^(S„ X [0, r]) n Vm{K) = ( U expi^([0, «] • {-MW)))) n Vra{K) 

(70) 

= T'^(Em,a) . 

Applying the coarea formula to the map ^, and using (77), yields 

(60) 



(60) r 

volT^(S„,a) < / ndet(Z)(H.,t)^)d(S^ xM) 

(Vl/,t)eS,„x[0,r] 

J Jo dety(>Vp) 



Changing the integration via the coarea formula applied to the map 11^ : M, and 

using (77), we get 

voir^(S^,a)< j j j\l + t^y''^ ■\c\iYiWp,t)\dtdYdp 

peM yeGr(TpM,m-l) ^ 

'■^ VolGr„_2,m-l I , ,M , , , 



where the expectation is with respect to Y chosen uniformly at random in Gi {TpM,m — 1). 
Using Theorem 4.2 with k = n — 2 and £ = m — 1, we obtain 



VOir-(S^,a)< J • E m-<^n-2-,m.e-'-^dtdp 



^'^-^ peM 

n-2 

^'^U''^ VOlGr„_2,„.-l • J2 • ^n,i{0^) ■ ■ ' VA^^) ■ 

i,j=0 
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So finally we get for rvol 7^(E^, a) = ^^J^^, using the formula (71), 

rvol T^i^m, a) < ^""[^"Z-''""' ■ J2 • • ■ ' Vj{K) 

VOi ^l-n rn 



i,j=0 



(71) 2m(n - m) fn/2\ \^ . j^. Oj ■ On-2-j v^immi r . x 

«' . ( V2) . g ■ ■ ° 

The remainder of this section is devoted to the proof of Theorem 6.1. In order to get a 
clear picture about where the convexity assumption M = dK with K G /C™(5"~^) comes 
into play, for the time being we will only assume that M C S""^^ is a compact oricntable 
hypersurface with a unit normal vector field u: M ^ S^~^. (Take M for example to be the 
smooth boundary of a non-necessarily convex domain of S'^~^.) 

We shall define a Sigma set UmiM) C Gr n,m, which extends the notion Sm(C) in the case 
M = dK with K G /C™(S'"^~"^) and C = conc(i^). As the general framework of Section 5.2 
and Section 5.3 suggests, it is useful to also define a lifted Sigma set S(Af) C 0(n) such that 
Eot(M) = n(S(M)), the canonical projection 11 : 0(n) GiCn,m being defined as in (65). We 
first define the lifted Sigma set S(M) before we give the definition of Em(M). 

In the following wc identify the tangent space TpM for p £ M with the linear subspace 
TpM = rwip)^ of M", which has codimension 2. If (p, Q) G F{M), then Q = [qi, . . . , qn-2) 
is an orthonormal basis of TpM and in the following we shall interpret Q as the matrix 
in M"^("-2) ^]jg columns Qi G M". Note that {p, Q, iy{p)) = (p, qi,..., qn-2, '^{p)) G 0{n). 
This defines the map 

$ : F{M) ^ 0(n) , {p,Q)^Q = {p, Q, u{p)) . 

We define the lifted Sigma set S(M) as the image of the map $ 

t ■= t{M) ■= ^{F{M)) = {g G 0{n) \ 3peM :Q-ei=p, Q-en = u{p)} , 

where G M" denote the canonical basis vectors. Analogous to the definition of we define 
for 1 < m < n — 1 the map (cf. (75)) 

: Gr(Af, m - 1) ^ Gr„,„ , {p,Y)^Rp + Y , 

where Y G Gr(rpM, m — 1) is interpreted as an (m — l)-dimensional subspace of W^. Further- 
more, we define 

:= E^(M) := $^(Gr(M, m - 1)) = {W^ G Gr„,^ \ 3p e M : p e W , W C u{p)^} . 

Note that we have n(I](M)) = S^, where 11: 0{n) GTn,m denotes the canonical pro- 
jection (cf. (65)), but we have a strict inclusion S(M) C n-i(E„(M)). Wc will see below 
(cf. Corollary 6.3) that S is in fact a smooth submanifold of 0(n) without posing any further 
assumption on M. In Corollary 6.5 we will see that S^(M) is a smooth hypersurface in Gr„^„ 
ifM = dK for K G ]C"^{S''-^). 
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The following commutative diagram provides an overview over the relations, which are 
central for the understanding of S^: 



F := F{M) — 
Hi, 

G :=Gr(M,m- 1) 



0(n) 



n 



(78) 



Gr 



where is defined as in (72), and is defined as the restriction of 11 (cf. (65)) to S. Note 
that if {p, Q) G F(M) with Q = (-71, • • • , 9n-2), then l>(p, Q) = Q = (p, gi, . . . , g„_2, z/(p)) € S 
and nb(p, Q) = (p, 1"), where F is the span of gi, . . . , qm-i- 

The tangent spaces of the fiber bundles F{M) and Gr(M, m — 1) have natural decompo- 
sitions into vertical and horizontal components (cf. Section 5.3) 



- (p,<9) 



(79) 



We denote the images of the vertical and the horizontal space of the frame bundle at (p, Q) 
under the derivative by 



(80) 



We will give next concrete descriptions of these spaces. 

Before wc state the next result about the derivative of $ we define the following linear 
maps for Q = {p, Q, z^(p)) G S 



E'" : Skew„_2 Skew^ 



en— 2 



Skewr; 



(81) 





( ° 





\ 




( ' 


T 

—a 





E'"{U) := 





u 







a 





-b 




1 





) 




\ 


6'^' 






where h := ■ Wp{Qa) . 



Here Wp : TpM — )■ 7],M is the Weingarten map of M at j5, cf. (8), and we recall that the image 
of € M"^("-2) equals TpM, as {p,Q,v{p)) e S. Note that E'" is an isometric embedding 
independent of the manifold M and the element Q & T,, while & is a linear injection, which 

depends both on M and on Q. For the sake of simplicity, we do not reflect this dependence in 
the notation. Note hirthcr that the images of and arc orthogonal subspaccs of Skew„. 

Proposition 6.2. Let ip,Q) G F{M) and Q := ^ip,Q) = {p,Q,v{p)) G S. Then we have 

T^t = Q ■ im{E'') , T^t = Q ■ im{E'') . (82) 

Proof. Since p = Qe\ and v{p) = Qcn, the image of the fiber {p} x F{TpM) under the map $ 
is given by 







( 1 








$(M X F{TpM)) = < 


Q- 





Q 









I 





1 



Q G 0(n-2) 
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This implies that the image of the vertical space T^^ q^F{M) under the derivative is given 
by T^t = Q ■ im(£:^), cf. (74). 

As for the horizontal space, let ^ G q^F{M) be represented by the curve cj? : M — )■ F(M), 
Ci?(0) = {p,Q) and cf{0) = ^, which is given in the following way: cp^t) = {c{t),Qt), 



I.e. 



where Qt is the parallel transport of Q along the curve c (cf. Section 5.3). Let us denote 
by C := c(0) G TpM the tangent vector at p defined by c. The image of the curve cp 
under $ is thus given by $(cir(t)) = (c(t), Qt, u{c{t))) =: Q{t), and the image of ^ is given by 
D$(^) = Q{0). It is sufficient to show that Q{0) is given by 





( " 




\ 




Q(0) = Q ■ 


a 





-b 


, where a = Q^C , ^' = Q^H^p(C) 




\ 




/ 





(83) 



To prove (83), note first that as Q{t) G 0{n) and (5(0) = Q, wc have (5(0) = Q ■ U with 
?7 G Skew„. Recall that the columns of Q form an orthonormal basis of TpM = p-^ Ci ^[p)-^, 



hence Q^TpM 



^n- Therefore, c(0) = C = Q<^ ^or some a G 



pn— 2 



which implies 



a = Q^C- The first column of U is thus given by C/ei = Q^Q{0)ei = Q^C, = By 
skew-symmetry of C/, this also gives us the first row of U. The zero matrix in the middle 
follows from the fact that the frame Q is parallel transported along c (cf. [8, §VII.3]). Finally, 
the last column of (5(0) is given by DpV^Q. This implies that the last column of U is given by 

Q^Dpi^ic) = -Q^mc) ■ 

As >Vp(C) G TpM = p-^ n vip)-^, the last column of U has the form (^-^^ ■ Therefore 

^p(C) = Q i^b^ = Qb, which implies b = Q'^Wp{(). The last row follows again by skew- 
symmetry of f/. □ 

Corollary 6.3. The map $ is an injective immersion of F{M) into 0{n). In particular, S is 
a smooth submanifold of 0(n) of codimension n — 1, and the tangent space TqT, at Q E 0{n) 
has the orthogonal decomposition TqT, = TqT © TqT. 

Proof. The fact that $ is a smooth injective map is obvious. Furthermore, by Proposition 6.2, 
it follows that the derivative at {p, Q) has full rank. The map $ is thus an injective 
immersion, and as the domain F{M) is compact, it is also an embedding. As for the dimension, 
we compute 

dim t = dim F(M) = dim M + dim 0(n - 2) = n - 2 + ("-^K"-^) = {n-2)(n-i) 
= dimO{n) — (n — 1) . 

The decomposition of the tangent space into the direct sum TqT = TqT® TqT follows from 

the decomposition of the tangent space of the fiber bundle (79). The fact that TqT and TqT 
are orthogonal follows from the description given in Proposition 6.2. □ 
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Now that we have a clear description of the hfted Sigma set S, we will transfer this 
description to via the projection map (cf. (78)), resp. its derivative DH^. Analogous 
to the lifted case we denote the images of the vertical and the horizontal space, cf. (79), of the 
Grassmann bundle at {p, Y) under the derivative by 



(84) 



Recall that the subspacc Skew^ defined in (66) serves as a model for the tangent spaces 
of Gr„^^. We consider the linear map 



where xi G 



a: 



D m— 1 



pn— 2 



Skewn 






T 
X2 








-Xl 












T 










\ 



(85) 



and X2 G 



•bn—m—l 



, and we denote its image by 



Note that Sn has dimension n 

(n— m— l)x (m— 1) _ 



Sn '■= im(a) C Skcw„ . 

■ 2. We now define the following linear maps for Q 



Skewr 



E' 



on— 2 



Sn 



(86) 









\ 










X 













/ 



E''{a) := 








T 

-ai. 










-bi 


02 








V 





where again b := ■ yVp{Qa), and a 



with ai, 61 G 



um— 1 



02,62 G 



pn— m— 1 



The map E^ is an isometry if the scalar product on Skew„ induced from M"^" is scaled by ^. 
Note that E"" is independent of M and Q. On the other hand, as in the lifted case, the map E^ 
depends on both M and Q. Unlike in the lifted case, the map E'^ need not be injective. We 
will see below in Proposition 6.4 that it is injective iff the restriction Wp,y of the Weingarten 
map Wp to the subspace Y = lin{(7i, . . . , qm-i} Q TpM, where Q = {p,qi, . . . , qn-2, ^ip)), has 
full rank. If this is the case, then im£''* = Sn- 

Note that we have an orthogonal decomposition of Skew„ into 

Skewn = im^;" e S'„ e MA^ , (87) 

where N G Skew„ was defined in (62). If we define l: m("-"^-i)x Skewn-2, via 
X i-> ( ^ ) , and if we denote by tt : Skew„ Skew„ the orthogonal projection, then we 
have the crucial relation 

E^iX) = E^tiX)) , E^{a) = n{E^{a)) . (88) 

Proposition 6.4. 1. Let {p, Y) G Gr(M, m-1) and W := ^m{p, Y), and letQ et^ be such 
that I[s{Q) = W. Then we have, using the model of tangent spaces of Grn,m described 
in Section 5.2, 



T^^m = {[Q,U]\U eim{E'')} , 



T^^m = {[Q,U]\U eimiEf')} 



(89) 
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2. The rank of the map is given by 

vkiE'') = n-m-l + vkWpx , 

where Wp^y : Y Y is the restriction of the Weingarten map Wp to the subspace Y 
ofTpM. Moreover, ifWpy has full rank and if we consider Sn to be endowed with the 
scalar product {Ui, U2) '■= \ ■ tr(C/^ • U2), then 

|det(£;'*)| = |dety(Wp)| , (90) 

where dety(>Vp) denotes the determinant ofWp^y (cf- Definition 4-i)- 

3. The rank of the derivative D(^p Y)^m is given by 

rkD(p y)$m = ^{n — m) — m + rkWp^y , 

Proof. (1) Since the derivative of the bundle projection Hj, (cf. (72)) maps the vertical spaces 
of the frame bundle to the vertical spaces of the Grassmann bundle, cf. Section 5.3, we have 

r^S„ D<^>^ (^fey)G') =^ Dn^mt) DU4Q ■ im(^^)) . 
Using DqIIs{QU) = [Q,U] (cf. (67)), we may continue as 

r^S^ = {[Q,U]\Ue 7r{im{E'))} = {[Q,U]\Ue imiE"^)} . 

Similar arguments show the description (89) of the horizontal space. 

(2) By the above choice of scalar product in 5„, the map a defined in (85) is an isometry 
on its image. Hence, instead of considering E^ we may focus on the map 

V«2 / \a2 

where as usual b = ■ yVp{Qa), and a = (^02)' ^ ~ (ft^) ^'^^^ f^ij^i G M™^^, 0-25^2 £ 
j^n-m-i_ j^g^ ^ g ]^(?i-2)x(ra-2) representation matrix of Wp with respect to the or- 

thonormal basis qi,. ■ ■ ,qn-2 of TpM, i.e., Aa = • yVp{Qa) = b. Furthermore, let A be 

decomposed as A = (^^ ^i, where Ai G and the other blocks accordingly. 

\A3 A4 J 

Then the representation matrix of a~-^ o is given by ( ). As Ai is the repre- 

V U ^n-m-lj 

sentation matrix of Wp,y with respect to the orthonormal basis gi, . . . , Qm-i of Y, we obtain 

rk(^'') = rk(a-i oE^) = n-m-l + rk(Ai) =n-m-l + rk(>Vp,y) . 
Furthermore, if Wp,y has full rank rk Wp^y = m — 1, then 

I det(£;'')| = I det(a-i o E^)\ = \ det(Ai)| = |dety(>Vp)| . 

(3) From (1) and (2) we obtain 

dimT^Sm = rk(S^) = (m - l)(n - m - 1) , 
dimT^E^ = rk(£;'*) = n-m-l + rk Wp,y . 
This implies that the rank of the derivative of is given by 

rkD(p,y)$^ = dimT^S^ + dimT^S^ = m(n - m - 1) + rk>Vp,y . □ 
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Corollary 6.5. Let M = OK with K G /C='"(5'^-i). 

1. Then is an injective immersion of Gi(M,m — 1) into G'rn,m, o,nd is a smooth 
hypersurface of Gr„^m. Furthermore, for W £ E^, the tangent space of at W has 
the orthogonal decomposition Tw^m = ^^^m © T^^^m- 

2. If W = Yij:{Q) G for Q G then the unit normal vector of E^ at W , 
which points into the component V^, is given by v-siW) = [Q,N] G TwGin,m, where 
N G Skew n was defined in (62). 

3. Let Um -.T^rn^M he given by W ^ p, where W r\ K = {p] (cf (76) ). Then the 
normal Jacobian of Um at W £ T,m is given by ndet{DwIlM) = dety(Wp)~^, where 
Y :=Wnp-^ e Gr(rpM,m- 1). 

Proof. (1) The injectivity of follows from Proposition 2.3. Furthermore, the Weingarten 
map Wp of M is positive definite for every p G M. Therefore, every restriction Wpx of 
Wp to a subspace Y G Gr(TpM, m — 1) is positive definite and in particular has full rank. 
It follows from Proposition 6.4 that is an injective immersion. By compactness of the 
domain Gr(M, m — 1) it follows that is an embedding. The image E^ is thus a smooth 
submanifold of dimension dim E^ = dim Gr(M, m—1) = m{n—m)—l. The claimed orthogonal 
decomposition of Ty^Em is a consequence of (89) and (87). 

(2) It is clear from the orthogonal decompositions Tw'^m = T^r'^m®T^^Tim and (87) that 
v-siW) lies in the orthogonal complement of T^i/E^ in Tjy Gr„^^. To see that it points into 
the component Vm, we consider the exponential map in this direction. We have 

exp^^(p • u^{W)) n(expQ(QiV)) U{Q • Qp) . 

The first column of Q ■ Qp is given by := cos p ■ p + sin p ■ u{p). As the normal vector u{p) 
points inside the cap K, the vector pp lies in the interior of K for p > small enough. In 
particular, exp^y(p • i^e(VF)) G Vm for p > small enough. 

(3) We may lift the function IIm : '^m — >■ M to the function IIm : E M, Q p, where 
Q = {p,qi, . . . , qn-2,J^{p))- So we have the following commutative diagram 




M P 



where W = lm{p, qi, . . . , qm-i}- We next show that the orthogonal complement of ker DwUm 
in Tw^m is given by the horizontal space T^^Y^ra- 

As the map is just the projection onto the first column, the kernel of -DqII M IS given 

by 

kcri^Qn,, = {[/ G TqE I QUei = 0} . 

It follows that TqYi C ]lct DqIIm and TgE H kcrDQllA.f = 0. Using the commutative di- 
agram (91), we obtain Ty^Y^rn ^ kerD^yn^ and Tij^^Urn H \ieT Dw^m = 0. The fact that 
Ty/Ym = T^^T^m © Ty^Yrn is an orthogonal decomposition implies that Ty^Ym = kev DwUm 
and hence T^E^ = (ker Dv^nM)"*"- 
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In order to compute Dw^m on the horizontal space T^S^, we consider for fixed Q = 
{p,qi, . . . , qn-2, ^{p)) € £ the following diagram 

B 

Mn—2 ^ Q ^ . rph V 
* >~>n * J-w^'m 

II DwYIm ■> 

M"-2 1 , 2],M 

where /3(J7) = [Q, U] for U e Sn, and 7(a) = YJi=i ■ Qi for a = (ai, . . . , a„-2f G M"-^. Let 
us check that this diagram is commutative: using 

DgU^iQ . E\a)) [Q,T,{E\a))] [Q,E''ia)] , 

we get 

DwnM,E\a)]) = DwIlj^iDQU^iQ ■ E^'ia))) DqU^{Q ■ E\a)) = 7(a) . 
Since /3 and 7 are isometric we obtain 

ndet(L>vKnM) = |det(S'^)|-^ |dety(Wp)|-i = dety(Wp)-^ , 

where the last equality follows from the positive definiteness of Wp. □ 

For the proof of Theorem 6.1 we make yet another definition in the lifted setting. For 
Q G E we define the direction i>{Q) G TQO{n) and the map S x M ^ 0(n) via 

u{Q) ■.= Q-N , ^{Q, t) = expQ(arctan(-i) • z>(Q)) , 

where N G Skew„ is, as usual, defined in (62). Abbreviating p := arctan(— t), we obtain 
from (63) that the map ^ is given by ^{Q,t) = Q ■ Qp, where Qp was defined in (63). 
Furthermore, denoting 11 := U.^ x id]R, we get the following commutative diagram 



S X M > 0(n) 



n n ■ (92) 

Em X M » Gln,m 

Proof of Theorem 6. 1 . The first part of the theorem, as well as the claim about the normal 
Jacobian of H^, is covered by Corollary 6.5. 

It remains to show the claim about the normal Jacobian of ^. We will compute the 

derivative via the lifting in the following way. Let W G and let Q = (p, Q, I'ip)) G 
E be a lifting of W, i.e., W = Uj:{Q). For (^, i) G T^w,t){^m x M) let i G Tq± be a lifting of ^, 
i.e., ^ = DIl^(^Q. Using the relations displayed in the commutative diagram (92), we get 

D^{^, i) = D^{DU{i, i)) DU{Dm i)) . (93) 
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Prom the explicit form "^{Q, t) = Q ■ Qp, where p = arctan(— t) and Qp defined as in (63), and 
^ = -(1 + tY\ and from = QpN, we get 



D^Q^t^^{o,i) = -{i + t')--Q-Qp-N 



Using (93), we thus have 



D^w,t)^io, 1) = L>n(-(i + t'r' -Q-Qp-N) 



(67) 



2\-l 



[Q-Qp,N\ 



(94) 



It remains to compute the derivative of ^! in the first component. 

Taking into account the orthogonal decomposition Tw'^m = ^^^m © T^r'Sm, cf. Corol- 
lary 6.5, we first consider the vertical space T^S^. By (89) every element ^ G T^^Hm is of the 
form ^ = [Q, E'"{X)] for some X G ■g^in-m-i)x{in-i) ^ Furthermore, a hfting ^ of ^ is given by 
= Q.^^(X), as DUj,{[Q,E''{X)]) = Q-E''{X), cf. (67). As for fixed p, the map Q ^ Q-Qp 
is linear, we obtain from ^'(Q, t) = Q ■ Qp 

DiQ,t)-^{i 0) = g • E^{X) ■ Qp Q-Qp- E^{X) , 

where the equality (*) follows from the fact that Qp only acts on the first and the last columns 
or rows. This implies via (93) and (67) 



0) = [Q-Qp, AE^iX))] = [Q-Qp, E-{X)] , 



(95) 



where, as usual, tt : Skew„ — > Skew^ denotes the orthogonal projection. 

As for the horizontal space, any element C, £ T^^T^m is of the form = [Q, E^{a)] for some 
a € M"-2, cf. (89). A lifting C of C is given by C = Q • ^''(a), as 

Z^n,(C) [Q,7r{E\a))] [Q,E\a)] = C • 
We thus get from ^{Q,t) = Q ■ Qp 

L>(Q,t)*(C, 0) = Q • E^^ia) ■Qp = QQp-Ql- E\a) ■ Qp . 
This implies via (93) and (67) 

D{w,t)^iC, 0) = [q • , TT (qJ • E\a) ■ Qp)] . (96) 
To finish the computation, note that we have for a,b ^ M""^ 





I 





T 

— a 







a 





-b 




\ 





6^ 








( 


—ca^ + sb'^ 


\ 




ca — sb 





—sa — cb 




K 


so'- + cb^ 


/ 



(97) 
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where we use the abbreviations s := sm{p), c := cos{p). Using the decompositions a = 02 ) ■ 

6 = (^^1 j with ai,6i G M"*-!, 02,62 G R'^'"^''^, we obtain 





/ 


( ° 


T 

—a 





TT 




a 





-6 






[ 










/ 





—ca2 + SO2 





(97) 







— sai — cbi 




ca2 — S62 










I 


saj + 



(98) 



Combining (98) with the formula in (96), we get an exphcit formula for the derivative of ^' 
on the horizontal space. 

To summarize the results from (94), (95), and (96) let us define Wp := Ii{Q • Qp), and set 

T^^ := D^{0 X M) , T^^ := D^{T^J^rn x 0) , T^^ := £'*(r/5,S^ x 0) . 

Recall the orthogonal decomposition Skew^ = MA/"© imE'" © im£''*, cf. (87). This implies the 
orthogonal decomposition 

Furthermore, the corresponding restrictions of D'^ yield a (reflection and a) dilation by (1 + 
t^)-! between x M and T^^ (cf. (94)), an isometry between T^S™ x and T^^ (cf. (95)), 
and a nontrivial linear map between Ty^T^m x and (cf. (96), (98)). This implies 



ndet(£>(vi/,t)*) = (1 + i^)"^ • I det L>*| . (99) 

We determine now | detD*|. 

We consider the following commutative diagram defining the linear map /x: — >■ M"'~^: 



TyyT,^ X 



pn— 2 



Wp 
P o a 

(m— l)+(n— m— 1) 



where ^: Skewn T^Sm, U ^ [Q,U], and a: M"-^ Skew^ as defined in (85). As a,/3 
are isometrics, we obtain 

|det:D^| = |det(/?oE'^)-i| • |det/i| = Idet^'*!"^ • |det^| dcty(Wp)-^ • |det/i| . (100) 
As for the linear map fi, wc obtain from (98) that n has the transformation matrix 
("^' + f-^ , ,),^('A = ( ^1 • M and 

V -sAg cIn-m-1 - SA4 J ' \b2 J V^3 ^02/ 

/ cAi + s/^-i CA2 \ / ai\ _ /sai + cbi 

\ -sAa c/n-m-i-sA4y ^02/ \ca2 - sb2 
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Using that | = tan p = —t and c = cos p = 1/ Vl + i^, we obtain 
det/x = detf^^^ + f-^ r '"^^ 

\^ -SA3 Cin-m-l - shi J 



t-^S In-m-1 + t^A 

= (l + t2)-("-2)/2.chy(Wp,t), 

where the last equality is a consequence of (51). From (99) and (100) we thus get 

r,Hptrn VI/^-n^/2^-l Idet/i| _ , 2wn/2 |chy(Wp,t)| 

ndet(Z^(H.,)*) - (1 + • - (1 + / • detyW, ' 

which finishes the proof of Theorem 6.1. □ 
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